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L Introduction 

1.1 Diagnosis of Breast Cancer 

X-ray mammography is the current gold standard in breast imaging and, in combination 
with a physical exam, detects breast cancer with a sensitivity of 85-95%. [1,2] Diagnosis of 
detected focal abnormalities is based primarily on morphological features such as boundary 
appearance, presence or absence of spiculations, and distribution and size of calcifications. Some 
lesions can be characterized specifically as either benign or malignant on the basis of these features 
alone; however, many are indeterminate, and require further evaluation. [3,4] Recent trials have 
shown that 70-80% of mammographically indeterminate lesions that progress to surgical biopsy are 
benign. [4-7] For women with particularly dense breast tissue and/or fibrocystic changes, the false 
positive rate for x-ray mammography can be as high as 90%. [5,8] For these women, the surgical 
biopsy serves only to provide a diagnosis since it is not necessary to remove a benign lesion unless 
it is causing pain. 

Recently, less invasive biopsy techniques have been developed which use a large diameter 
needle fitted with a cutting device to extract cores of tissue from a lesion for histological analysis. 
In the case of a solid, palpable lesion, core biopsy is sometimes used to provide a definitive 
diagnosis even in the case of a negative result; however, this application is still very controversial. 
[9-12] Although the reported accuracy of this core biopsy technique in small studies has improved 
rapidly since its inception (largely because of improved sampling accuracy using needle guidance 
by imaging methods), a conservative approach demands more extensive testing in large populations 
with surgical confirmation before a meaningful evaluation of its reliability can take place. [13] 

Development of accurate, non-invasive methods for characterizing mammographically 
indeterminate lesions represents an active area of current research in medical imaging. Ideally, a 
new diagnostic test would have reliably low false negative rates and would significantly reduce the 
number of surgical biopsies for benign lesions currently being performed; however, even if the 
false-negative rate associated with the test were unacceptably high, a non-invasive test with high 
positive predictive value for malignancy would have a profound impact on current surgical practice 
for breast cancer. At present, patients undergo surgical biopsy, then return for a second surgical 
procedure if the biopsy is positive, allowing time for patients to assess their treatment options in the 
intervening weeks. A non-invasive predictor of malignancy could eliminate this two-step 
procedure, reducing the risk of disease progression and the psychological stress of delaying 
treatment. 

1.2 Patient Prognosis 

In 1990, Fearon and Vogelstein published their work elucidating the genetic steps to 
malignancy in colon cancer. [14] It is now widely believed that breast cancer also evolves in a 
similar way from normal ductal epithelium, to hyperplasia, to atypical hyperplasia, to carcinoma in 
situ and finally to invasive carcinoma. [15] At each stage in this progression away from normal 
epithelium toward malignancy and invasion, the prognosis for the patient worsens. Following a 
surgical biopsy, the extent of malignant transformation is determined and the metastatic potential of 
a tumour evaluated to best determine clinical treatment. 

The recent technical improvements in mammography and increase in the number of women 
who undergo screening have dramatically shifted the stage of cancers seen at pathology. 
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Approximately 66% of patients with newly detected breast cancer will have no detected lymph 
node involvement, and 20-30% of these newly-detected lesions will be found while still confined 
within a duct. [4,6] Fifteen to twenty per cent of node-negative patients with invasive carcinoma 
will relapse in the 10 years following surgery. [16] Currently, there is no proven technique for 
identifying this high-risk subpopulation within the node-negative group. For these women, 
standard pathological indicators such as tumour size, histologic differentiation, and hormone 
receptor status have shown variable correlation with disease recurrence and patient survival. [17] 
For these localized breast cancers, surgical removal of the tumour supplemented by radiation 
therapy has replaced mastectomy for many women. [18,19] The development of prognostic 
indicators to identify patients from this group who are likely to benefit from systemic 
chemotherapy is also an active and important area of current research. Currently, there are no 
accepted guidelines for clinical treatment of women with ductal carcinoma in situ (DCIS). [20] 
Clinical care of patients with early-stage breast cancer is controversial, and the development of a 
test to assess the metastatic potential of these lesions is an important area of research in breast 
pathology. 

A large body of recent work exists concerning the development of new prognostic 
indicators for early-stage carcinomas. The list of new markers currently being evaluated is long 
and growing. New markers which show promise are those reflecting cellular characteristics such as 
DNA content, ploidy (define), cell cycle status, and proliferation indices. In addition, new genetic 
tests are being created to evaluate the expression of oncogenes. For example, the tumour 
suppressor p53 gene has been receiving much attention as a prognostic indicator. Overexpression 
of the p53 mutant form has been found to correlate with high proliferative activity and disease 
recurrence in some studies. [21] Unfortunately, these relationships have not been found in all cases, 
[22] and this new area of research is in need of more concrete results before these indicators can be 
used routinely for clinical decision-making. None of the multitude of currently available prognostic 
tests are sufficiently reliable for predicting relapse or survival, and decisions about patient care 
must be made by weighing many prognostic indicators, in addition to assessing individual patients’ 
needs. 

1.3 Tumour Angiogenesis 

Tumour progression is known to be dependent on the ability to stimulate the growth of new 
blood vessels. [23] This process, known as tumour angiogenesis, is integral to all solid malignant 
neoplasms, and is required in order to supply nutrients and oxygen to the expanding tumour. In 
tumour spheroids, transport of nutrients and oxygen by diffusion was shown to limit tumour growth 
to a size of 1-2 mm diameter. [24-26] In a chick embryo model, however, implanted tumour cells 
were shown to stimulate the growth of new blood vessels, and to thereby achieve rapid growth 
beyond this diffusion-limited size constraint. [27] It was therefore believed initially that the switch 
to an angiogenic phenotype would occur when a tumour reaches its maximum diffusion-limited 
growth potential. More recent results have indicated that the onset of angiogenesis can occur much 
earlier in tumour progression, and may even occur well before the switch to a malignant phenotype. 
[28] Recent studies evaluating microvessel density around in situ breast carcinoma have shown 
increased vascularity even around these very early stage lesions. [29-31] More strikingly, cells 
extracted from premalignant lesions have demonstrated the ability to induce new vessel formation 
in model tissue systems. Using a rabbit cornea assay, cells extracted from atypical hyperplasia and 
even from histologically normal lobules from cancerous breasts have been shown to stimulate 
angiogenesis. [32-35] Although these results are preliminary, and mechanisms of angiogenesis are 
still not fully understood, it is evident that even very early malignancies, and some premalignant 
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lesions undergo marked changes in vasculature. 

1.3.1 Implications for Diagnosis : It has been hypothesized that the characteristic new vessel 
growth associated with tumour angiogenesis could allow differentiation between malignancies and 
benign breast lesions using medical imaging techniques. [36-40] Specifically, it has been speculated 
that the number of small vessels per unit volume of tissue (microvessel density) would be higher 
for malignant lesions. Some imaging researchers have even asserted that benign lesions are not 
expected to be vascularized at all, and that the detection of blood flow around a lesion would 
therefore strongly indicate malignancy. Although many attempts premised on this assumption have 
been made to diagnose breast disease using IR thermography, positron emission tomography, 
Doppler ultrasound, and magnetic resonance imaging techniques, none of these studies has 
achieved the highly reliable differentiation between benign and malignant lesions that was 
expected. This has led some authors to suggest that the absolute number of microvessels in benign 
and malignant lesions may overlap in some cases. More hopefully, several recent imaging studies 
have suggested that an increased number of these vessels specifically in the lesion periphery may be 
a characteristic that is uniquely associated with malignancy. [41,42] 

To date, most studies evaluating micro vessel density in breast using histological techniques 
have focussed on invasive breast carcinomas, and the increased number of blood vessels associated 
with a strong tendency to metastasize. Although comparisons of microvessel density to normal 
breast tissue have been made, relatively few of these studies have included benign breast lesions, 
and very little data about vascularity in benign tumours is actually available. The hypothesis of 
benign lesions being avascular appears to have been based almost solely on what is known about 
breast cancers, extrapolating through “less invasive” to “benign”. If the number of blood vessels or 
their spatial distribution around malignant lesions is, in fact, characteristic for malignancy, a non- 
invasive measure of blood flow around ein indeterminate lesion would add valuable information to a 
diagnostic work-up; however, histological studies evaluating distributions of vessels in benign 
lesions are certainly required to test this hypothesis. 

1.3.2 Implications for Prognosis : The entry of cancer cells into the bloodstream and subsequent 
formation of distant metastases is facilitated by the increasing number of small new blood vessels 
around the tumour as it grows. Studies of patients with invasive breast carcinoma have shown that 
the number of microvessels in angiogenic “hotspots” in the tumour periphery correlates strongly 
with the extent of metastatic progression. [43-47] In one study, each 10-fold increase in the 
number of microvessels per x200 microscope field corresponded to a 1.6-fold increase in the risk of 
metastasis. [43] The number of microvessels in these hotspots has also been shown to outperform 
l)nnph-node positivity as a predictor for the development of distant metastatic disease. [44-46] 
Moreover, micro-vessel density is the only prognostic indicator that is a statistically significant 
predictor of overall survival for node-negative women. [44,47] A non-invasive diagnostic 
technique that measures microvessel density in invasive carcinomas would identify node-negative 
lesions with high metastatic potential and would immediately indicate which patients are most 
likely to benefit fi-om adjuvant systemic therapy. 

A measurement of microvessel density may also have value in prognosis for in situ 
carcinoma and for premalignant disease. It has been proposed that the degree of angiogenic activity 
around premalignant lesions and early-stage cancers may be an indicator of which lesions will 
progress to more invasive disease. [48] Preliminary results fi'om a retrospective study by 
Guinebretiere [49] showed that vessel density around hyperplastic lesions correlated with risk for 
subsequent development of invasive carcinoma and, further, that vessel density had a higher 
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prognostic value than atypia. A non-invasive method for measuring microvessel density around 
premalignant lesions and in situ cancers may provide early identification of lesions that have the 
potential to progress to more invasive disease, and require more aggressive therapy, or long-term 
patient follow-up. 

1.4 Magnetic Resonance Imaging (MRIl 

Breast imaging using MRI has several advantages over x-ray mammography: (a) image 
contrast can be based on one or more of several independent tissue parameters (Tl, T2 and water 
content) and can therefore be manipulated to optimize contrast between lesions and normal 
parenchymal tissue, (b) optimal display of tissue close to the chest wall can be obtained, since no 
compression of the breast is required, (c) images are obtained either as 2-dimensional slices in any 
orientation, or as 3-dimensional volumes, eliminating overlap of internal breast structures on the 
image. 

The focus of research in MR imaging of the breast has been on contrast-enhanced studies. 
It was originally believed that these contrast exams would distinguish malignant lesions from 
benign tumours with a high degree of reliability, since microvessels associated with malignant 
lesions were thought to be more numerous and more permeable than in normal tissue. Initially, 
simple snapshot “before” and “after” images were acquired, and areas of focal enhancement were 
evaluated for percentage signal increase. These MR studies reported detection of breast cancer with 
high sensitivity; however, diagnostic specificity was widely variable across different studies (20- 
85%). [50-53] Most solid malignancies showed marked contrast enhancement, allowing easy 
differentiation from mammographically similar features such as cysts, old scar tissue, and dense 
normal breast tissue. Unfortunately, some benign diseases such as atypical hyperplasia and 
fibroadenomas showed similar absolute enhancement to malignancies, and some malignancies 
demonstrated benign-type enhancement. 

Despite widespread efforts to refine MR techniques and analysis of these contiust-enhanced 
studies, little improvement in specificity resulted. [39,42,54-56] Extremely high sensitivity to 
cancer using this method is reported (94-100%); unfortunately, a large number of false positives 
has also been reported, where the majority of enhancing benign lesions are fibroadenomas or 
fibrocystic change. 

Recently, new MR hardware has become commercially available that allows higher spatial 
resolution imaging (300-700^m in-plane resolution with 2-3 mm slice thickness). This improved 
resolution allows more detailed evaluation of lesion borders, and internal structure. In an initial 
study of patients with a palpable or mammographically visible mass, a high spatial resolution 
image was obtained before contrast agent injection, a series of lower spatial resolution images was 
acquired rapidly during the first 4-7 minutes post-injection to obtain high temporal resolution 
enhancement curves, and a second high spatial resolution image was then acquired using identical 
parameters to the first image. [41] These preliminary high spatial resolution contrast imaging 
results have indicated that the additional information about architectural features combined with 
initial slope calculations may result in increased specificity for breast cancer. Specifically, it was 
found that malignancies tended to demonstrate higher contrast enhancement around the lesion 
periphery and also tended to have spiculated borders, whereas benign lesions such as 
fibroadenomas exhibited a more diffuse, uniform enhancement with smooth borders. In a further 
study looking specifically at patients with DCIS, this technique detected ten of thirteen DCIS 
lesions, where three of the ten identified lesions were not visible on mammography. [57] These 
results have demonstrated considerable promise for combining information about tumour blood 
flow with high-resolution imaging techniques. 
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Although initial high-spatial resolution contrast-enhanced imaging results are promising, a 
fundamental problem with contrast-enhanced methods is that the signal enhancement seen reflects 
two inseparable aspects of tumour blood vessels: total vessel surface area within a voxel, and 
vessel permeability. Moreover, contrast exams are expensive and contrast agents are not well- 
tolerated by all patients. A non-invasive technique that is sensitive to only one aspect of the 
vasculature e.g. blood flow rate provides simpler interpretation. 

The goal of this thesis was to test the feasibility of using the Intravoxel Incoherent Motion 
Imaging technique (I VIM) to measure microvessel blood flow around human breast lesions. I VIM 
is essentially a modification of an established protocol used for over thirty years to measure small 
water motions due to diffusion. In IVIM imaging, the same motion-sensitive protocols are used, 
where the protocols are tuned for sensitivity to the motion of the water component of blood flowing 
in small vessels. The IVIM technique has been demonstrated to be sensitive to blood flow in 
microvessels in the brain [58], and provides two quantitative measures: 1) absolute determination of 
the average blood flow velocity in these vessels within each tissue voxel, and 2) a relative measure 
of the total amount of blood flowing in these vessels at the time the image was acquired. It is 
possible that the measurement of these parameters in breast tissue may provide valuable 
information about microvessel densities and distributions around breast lesions. 

The thesis outline was as follows: 1) The hypothesis that spatial distributions of blood 
vessels around a solid breast lesion are specific for malignancy was tested in a study of archival 
biopsy tissue stained for blood vessels; 2) The IVIM method was evaluated using a mouse tumour 
model; 3) Hardware necessary for the application of the IVIM technique (or diffiision imaging) to 
human breast cancer in a clinical setting was designed, constructed and tested for feasibility. 

2. Distributions of Micro vessels in Breast Lesions 

A comparison study of spatial distributions of microvessels in human breast carcinomas and 
in fibroadenomas was carried out. Fibroadenomas were selected from other forms of benign breast 
disease for comparison to carcinomas because MR contrast enhancement patterns can often mimic 
those seen in invasive cancers. Moreover, fibroadenomas are the most commonly occuring benign 
solid lesion, and constitute 30-50% of benign lesions seen at surgical biopsy. [59] 

Cases were selected from the files of the Department of Pathology at the London Health 
Sciences Centre, and were stained using Factor Vlll-related antigen to identify endothelial cells. 
Microvessels were counted on x200 magnification in 10 fields in the lesion periphery, and 10 fields 
in the lesion centre. Where normal tissue was present, microvessels were also counted in 5 fields. 
The lesion periphery was easily identified on fibroadenomas as an enclosing fibrous capsule, and 
was defined as the leading edge of invading cells for the invasive carcinomas. Vessels were 
identified according to the criteria of Weidner et al. [43] where even a single stained endothelial cell 
was counted as a vessel. 

A number of interesting results were obtained from this work. Many of the fibroadenomas 
in our sample contained, on average, more microvessels than some of the invasive carcinomas we 
looked at. This finding directly contradicts the assertion by imaging researchers that benign breast 
lesions are not vascularized, and also explains, to some extent, the overlap in imaging 
characteristics seen for fibroadenomas and breast cancers. On a more optimistic note, our findings 
revealed that the numbers of small vessels were higher around the boundaries of invasive 
carcinomas, whereas these small vessels were more uniformly distributed in fibroadenomas. This 
characteristic spatial distribution of blood vessels may allow differentiation between fibroadenomas 
and carcinomas using non-invasive imaging methods sensitized to flow in microvessels. A 
secondary, but interesting observation was that for the three women aged forty or older who had 
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fibroadenomas, microvessel density in the surrounding normal breast tissue was extremely high 
compared with all other women in the fibroadenoma group (for the women with breast cancer, no 
“normal” tissue was present). 

A manuscript describing this work has been accented for publication in the journal 
Radiology as: K. Weind, C.F. Maier, B.K. Rutt, M.Moussa, “Comparison of Micro vessel 
Distributions in Invasive Carcinomas and Fibroadenomas of the Breast” (included as Appendix 1). 
Portions of this work have also been presented at the Annual Meeting of the United States and 
Canadian Academy of Pathology, Orlando, FL, March 1997; and at the 1997 Annual Meeting of 
the Radiological Society of North America (RSNA). 

3: Diffusion Imaging in Implanted Human Breast Tumours 

Quantitative diffusion measurements were performed in tumours arising from inoculation of 
nude mice with human breast cancer cells to determine values for the Apparent Diffusion 
Coefficient for various components of solid tumours. Diffusion maps were compared to histology, 
and correlated well with gross tumour morphology. These diffusion measurements were carried 
out to provide a foundation for the implementation of IVIM experiments in these tumours, by 
characterizing the behavior of the signal arising from the static spins within the tissue. While these 
experiments were undertaken as a stepping-stone to IVIM imaging, the results are significant for 
researchers in cancer therapy, and show that the degree of necrosis resulting from a particular 
therapy can be monitored non-invasively using MRI. 

These experiments were carried out under the supervision of Dr. Degani, Dr. Rutt, and Dr. 
Bendel. All experiments were performed at the Weizmann Institute of Science in Rehovot, Israel in 
the lab of Dr. Degani, and this work was completed prior to funding by the U.S. Army Predoctoral 
Fellowship award. In addition to providing details related to water diffusion in tumours necessary 
for planning the IVIM experiments, this work was highly valuable as a learning experience. By 
performing these experiments in Dr. Degani’s lab, I was able to benefit from her years of 
experience working with this tumour cell line in nude mice, and also from Dr. Bendel’s MR 
imaging expertise. 

Ideally, IVIM experiments would have also been undertaken in Dr. Degani’s lab; 
unfortunately, the limited length of my stay there (4 1/2 months) did not allow for the extensive 
software and hardware testing that is required for such sensitive work. Diffusion imaging alone is 
highly demanding on MR system performance, and the ability to reliably extract information about 
5-10% of the collected signal requires thorough and meticulous system evaluation. We chose 
instead to continue the work here in London, Ontario, on our newly acquired 4 Tesla whole-body 
imager. 

A paper describing this work Cincluded as Appendix 2) has been published as; C.F. Maier, 
Y. Paran, P. Bendel, B.K. Rutt, H. Degani, "Quantitative Diffusion Imaging in Implanted Human 
Breast Tumours". Magnetic Resonance in Medicine, 37, pp.576-581, 1997. This work was also 
presented at the 4th Meeting of the International Society for Magnetic Resonance in Medicine, New 
YorkCity,NY, April 1996. 

^ IVIM Imaging in Implanted Tumours 

In 1986, Le Bihan recognized that under certain conditions, blood flow in randomly 
oriented capillaries could be considered as a “pseudo-diffusion” process and proposed using the 
PGSE method to visualize this microcirculatory. [58] This application of the PGSE sequence to 
imaging flow in capillaries was called Intravoxel Incoherent Motion Imaging, because the condition 
for this flow to mimic a diffusion process is that the flow within a voxel be “incoherent” with 
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respect to the imaging time. 

For voxel-sizes that can be obtained clinically (lmm^-10mm^\ we expect to find thousands 
of capillaries, with a more or less random distribution of orientations. Flow in these vessels 
therefore resembles a random-walk or pseudo-dififusive motion which can be characterized by an 
ADC. The relationship of this ADC to the vascular geometry is governed by the flow velocity of 
water molecules inside the vasculature. If the length of time between diffusion-encoding pulses, A, 
is small, the average velocity can be derived directly from the ADC and is not dependent on the 
capillary network geometry. To ensure that this requirement is met, it is important to use high 
gradient strengths for large values of b (rather than increasing b by extending A). Additionally, the 
use of large, short g pulses allows imaging with short TE, so that little signal loss occurs from T2 
decay during the imaging sequence. 

The total signal arising from a voxel consists of a contribution from spins flowing in the 
microcirculation (5-10% of total signal) and from spins undergoing Brownian diffiision in the tissue 
(90-95% of total signal). This signal will exhibit a biexponential decay when plotted against 
increasing values of b. The fraction of all MR-visible spins within a voxel which are moving in the 
microcirculation can be determined experimentally from this curve by comparing the amount of 
signal attenuation due to fast decay (area under the exponential characterized by D*) to the amount 
of signal attenuation from the “static component” (area under the exponential characterized by D, 
the ADC describing tissue water diffusion). Although this D* decay is due to the spins that are 
moving within the microcirculation, and does not contain a contribution from static spins inside the 
vasculature, it is expected that the fraction of moving spins determined will be related to the density 
of capillaries within the tissue voxel. 

Since the publication of Le Bihan’s work, much controversy has surrounded the I VIM 
theory, and Le Bihan’s experimental results. The strongest evidence in support of I VIM has been 
the observation of biexponential signal decay curves, and the derivation of biologically reasonable 
flow fractions and velocities from these curves. [60-65] Unfortunately, these results have not been 
reproduced reliably in all studies undertaken to test this technique. [66,67] Others have pointed out 
the inherent difficulties of fitting data to a biexponential function, and estimated that the signal-to- 
noise in these experiments was not high enough to derive flow parameters reliably. [68] Many 
investigators have argued that the source of the biexponential decay is more likely to be 
macroscopic motions, such as involuntary patient movement, pulsatile brain motion, or CSF flow. 
[69,62] Others have taken issue with Le Bihan’s proposed relationship between the 
pseudodiffusion coefficient, D*, and the microvascular geometry and flow parameters, and have 
developed more complex models to describe the decay curves. [70,71] 

The goal of this study was to perform preliminary investigations evaluating the possibility 
of implementing PGSE-based IVIM for imaging breast tumours. In these experiments, we 
attempted to control as many of the potential problems with IVIM measurements as possible by 
using an animal model and MR system in which the vascular signal and signal-to-noise could be 
maximized, and macroscopic movements minimized. First, to maximize the size of the 
microvascular signal, we chose a tumour model that was known a priori to be highly vascularized. 
Secondly, mice were anaesthetized using an inhalation anaesthetic, isofluorane, that is known to 
have minimal effect on tumour blood flow. An rf coil was specially constructed to give excellent 
B, uniformity and signal-to-noise, and further, extensive system characterization was performed 
using phantoms to ensure optimal system performance. 

4.1 Methods 

4.1.1 Evaluation of System Performance 
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In normal tissue, the number of spins within the microcirculation contributing to an MR 
image is a small fraction (< 5%) of the total number of MR-visible spins. This fraction is slightly 
higher for well-vascularized tumours, but is not expected to exceed 10%. Because the total MR 
signal arises from a superposition of the signals from the spins that are flowing in the vascular 
compartment and the spins that are undergoing molecular diffusion in the tissue, any experiment 
attempting to isolate a signal from the “flowing fraction”, f, requires careful attention to 
experimental design details and rigorous system testing. 

One of the most significant causes of inaccuracy in diffusion imaging are the cross-term 
contributions from imaging gradients and diffusion gradients. [72,73] These terms were minimized 
for our pulse sequence by moving the readout compensation pulse from its normal location before 
the 180° pulse to immediately before the acquisition window. The size of the remaining cross¬ 
terms was evaluated in the following way: two sets of PGSE images of a water phantom were 
acquired for all possible image orientations, with five diffusion-weightings ranging between b = 1 
s/mm^ and b = 450 s/mm^, and each set acquired with both positive and negative diffusion 
gradients. The difference in D calculated with images acquired using positive gradients and with 
those acquired using negative gradients is a sensitive indicator of the presence of imaging gradient 
cross-terms. 

A standard PGSE sequence was used with TE/TR = 56ms/2000ms, A=30ms, 6=13ms (A 
and 6 are defined as in Figure 1). Images were acquired on a Varian/Siemens Unity INOVA 4 
Tesla whole-body MR imager equipped with 25 mT/m gradient coils (the maximum gradient 
strength available at the duty-cycle required for our pulse sequence was 22 mT/m), using a custom- 
built 4.5 cm diameter quadrature birdcage rf coil. For our PGSE experiments, gradients were 
applied simultaneously along all three axes, thereby increasing the diffusion gradient strength by a 

factor of VJ to 38 mT/m. The maximum available gradient strength, together with the longest A 
possible for this TE gave a diffusion weighting of b = 450 s/imn^ A 5 cm field of view was 
qualitatively judged to give adequate spatial resolution for depiction of anatomical detail in the 
mouse. The choice of a 128 x 256 matrix size (6 mm thickness) with TR = 2 s required 
approximately 4.5 min per diffusion-weighted image. For all experiments, the rf coil and mouse 
(or phantom) were mounted on a separate, distinct holder which was completely mechanically 
isolated from the rest of the system. The function of this holder was to isolate the volume being 
imaged from any system noise due to gradient coil or coldhead vibrations. [74] 

To test the stability of the system over an approximately one-hour scan time, preliminary 
tests with a doped-water phantom were performed. The water phantom used in these experiments 
consisted of a 20cc plastic syringe (approximately 2 cm in diameter) that was filled with Immol 
MnCl in 0.9% saline. Images of the phantom were acquired continuously for 1.25h (with the 
parameters described above), toggling back and forth between minimum and maximum diffusion- 
weightings: g = 1 mT/m (b = 1 s/mm^), and g = 38 mT/m (b = 450 s/mm^). The phantom was 
maintained at 29 +/- 0.5 °C throughout the experiment by the use of several electrical resistance 
heating pads which lined the main magnet bore at the location of the rf coil. 

Images were transferred to a Sun Microsystems ULTRA SPARC 1/140 workstation for 
analysis using MATLAB software Version 5.0 (The Mathworks, Inc. Natick, MA, USA). Mean 
signal intensities and standard deviations were measured inside a ROI defined in the phantom. 

4.1.2 MRl Measurements in a Tumour Model 

B16F1 murine melanoma cells were cultured as described previously. [75] Four- to six- 
week-old female nude mice were injected subcutaneously in the flank with 1x10^ cells in a volume 
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of 0.1 ml of medium. This cell line is known to form well-vascularized tumours rapidly. Two 
groups of six mice each were injected with cancer cells. Tumours were allowed to grow for two to 
three weeks after injection of cells to sizes of approximately 0.7 to 1 cm in diameter. Mice were 
provided with food and water ad libitum. 

Isofluorane, an inhalation anaesthetic known to have minimal effects on the cardiovascular 
system and tumour blood flow, [76] was used to anaesthetize the mice during imaging. To 
anaesthetize a mouse for set-up prior to the MRI exam, the mouse was enclosed for approximately 
10 to 15 min in a container filled with 1 to 2% isofluorane vapour in oxygen. The mouse was then 
placed on a small acrylic holder, and paper surgical tape was used to hold the mouse in position 
during the MR exam, with particular attention paid to restraining the tumour-bearing flank without 
impairing breathing, or blood flow to the tumour. The mouse and holder were then placed inside a 
4.5 cm diameter quadrature birdcage coil built on an acrylic cylinder. The ends of this cylinder 
were sealed, and a mixture of oxygen and isofluorane vapour was circulated through the coil. 
Normal thermia was maintained during the experiments using the resistance heating pad 
arrangement described above. After initial set-up and localization images, transverse-orientation 
PGSE images of a tumour slice were obtained for approximately 1.25h, according to the pulse 
sequence protocol described above in §4.1.1. 

4.2 Results 

4.2.1 Evaluation of System Performance 

Diffusion measurements obtained in the water phantom were extremely accurate and 
reproducible. Exponential diffusion decays were obtained using five diffusion-weighted images for 
each image orientation, and gradient polarity. These curves were plotted for each imaging slice 
fi-om the mean signal intensities within an ROI defined inside the phantom. Diffusion coefficients 
were calculated for all slices fi:om a least-mean-squares fit of the natural logarithm of these data to 
the diffusion b factors. Results are shown in Table 1 for calculated diffusion coefficients for all 
possible image orientations and both diffusion gradient polarities. 


Table 1 Measured Diffusion Coefficients in a H 2 O Phantom at Room Temperature 


Slice Orientation 

Diffusion Gradient Polarity 

D (xlO-' mmVs) 

transverse 


2.18+/-.01 

transverse 

- 

2.20 +/- .01 

transverse 

gm^ 

2.19+/-.01 

sagittal 

4- 

2.13 +/- .01 

sagittal 

- 

2.22 +/- .01 

sagittal 

gm 

2.18+/-.01 

coronal 

+ 

2.14 +/- .01 

coronal 

- 

2.22 +/- .01 

coronal 

gm 

2.18+/- .01 


^ D calculated from geometric mean of signal intensities according to the method of Neeman et 
al. [72] 
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In initial experiments using the doped-water phantom placed on the patient table, 
reproducible results could not be obtained. The signal intensity was inconsistent from scan to scan 
and variations as large as +/- 50% were observed in the measured D. The source of these errors 
was attributed to vibrations of the patient table, which is indirectly mechanically coupled to the 
gradient coils through the table supports. These effects were most readily observed in highly 
diffusion-weighted images (i.e. large gradient amplitudes) as ghosting in images and strong 
smearing in the phase-encode direction. This kind of structured system noise has been observed 
previously by Conturo et al. [74] on another Siemens system with a similar patient table 
configuration as the 4 Tesla scanner used in this study (Siemens Magnetom Vision 1.5 Tesla 
scanner). 

Figure 2 shows results from the 1.25 hour scanning protocol described in §4.1.1. In these 
experiments, the phantom was movmted in isolation from the system using the custom-built holder. 
The data plotted corresponds to the mean signal intensity inside a ROI defined over the whole 
phantom (error bars = +/- SEM). The small standard errors (which are on the order of the size of 
the plotting symbols), and temporal stability of these signal intensities indicate the extremely high 
level of system performance achieved with this scanner. The diffusion coefficient calculated from 
these measurements, 2.45 x 10'^ mmVs agrees with literature values for water at 29°C. [77] 

4.2.2 MRI Measurements in a Mouse Tumour Model 

Images of tumours in four different mice were acquired according to the protocol described 
in §4.1.2. Although typical signal-to-noise on these images was excellent, (S/N « 400:1 for b = 1 
s/mm^, compared to S/N « 250:1 for b = 450 s/mm^ where noise was estimated as the standard 
deviation in a ROI in the background), phase smearing from motion was obvious. The signal 
intensity measured in a region of obvious phase smearing ranged from approximately 5% to 25% of 
the signal intensity within the tumour for the low and high b factor images respectively. In addition 
to motion artifacts in the images, more significant motions over the course of the exam were 
evident from visual inspection of the images. Figure 3 shows a typical timecourse of measured 
signal intensities from a mouse tumour for the 1.25 hour scanning protocol. The data plotted 
correspond to the mean signal intensity inside a ROI defined within Ihe tumour (error bars = +/- 
SEM). Comparison to Figure 2 demonstrates the inaccuracy of the results obtained for the mouse. 
This variability is easily explained by the presence of significant motion ghosts and phase-smearing 
in the images. ADCs calculated from consecutive pairs of these images varied by approximately 
13%. 

4.3 Discussion 

PGSE imaging using a doped-water phantom produced artifact-free images with extremely 
high reproducibility and accuracy. In Table 1, calculated diffusion coefficients are given for all 
three image orientations. Because the diffusion gradients are applied along all three axes, the 
diffusion coefficients for positive or negative gradient polarities should be the same regardless of 
the image orientation. The small differences observed (less than 3%) are likely to be due to eddy 
currents. Eddy currents are different for the three gradient coils, and also depend on the direction of 
the applied current. They therefore contribute differently to the diffusion-weighting b-factor for 
each of the imaging orientations and gradient polarities through their interaction cross-terms with 
the imaging gradients. 

For mouse imaging, macroscopic motions during the image acquisitions resulted in 
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unacceptable phase-smemng and image-to-image misregistration. The mouse tumour model was 
chosen for our experiments because of the prior, published data from other groups studying the 
effects of isofluorane on blood flow in solid tumours in mice. The effect of anaesthetic can vary 
dramatically for different animals, and additional calibration experiments would have been 
necessary for a different tumour model (using an independent measure such as laser Doppler 
anemometry). Moreover, it is unlikely that a different choice of animal tumour for these 
experiments would provide significant improvements over our model, since small macroscopic 
motions will be a problem for any animal study using such highly motion-sensitive pulse 
sequences. In general, motion during an imaging sequence that requires more than one excitation 
for all the phase-encoding steps produces random phase shifts that can add or subtract to the effect 
of the phase-encoding gradient from line to line in k-space. These can result in severe ghosting 
artifacts and smearing of the signal intensity in the phase-encode direction. In addition, the use of 
strong gradients that increase sensitivity to small velocities means that even subtle motions during 
the scan are an important source of error and image artifacts. 

Because macroscopic motions during a scan can often not be avoided, several technical 
strategies have been developed for overcoming these problems. Two main approaches have been 
extensively tested and proven effective: motion-correction techniques, and rapid imaging 
techniques that can acquire an entire image in less than one second (e.g. echo-planar imaging). One 
motion-correction approach is to use first-order flow-compensation to zero the phase artifacts from 
linear motion; however, this requires a significant reduction in diffusion-weighting. [78] Another 
approach is the use of a second, non-phase-encoded “navigator echo” in the PGSE sequence to 
correct phase errors due to rigid linear and rotational motions. [78,79] The chief limitation of this 
technique is that it is unable to correct for higher orders of motion, or motion that is spatially 
varying across the slice. Moreover, it is unlikely that the addition of navigator echoes can fully 
correct the motions observed during our long exam times. In echo-planar imaging (EPI), the entire 
set of echoes necessary to form an image can be collected with only a single if excitation pulse 
(single-shot EPI). The complete image acquisition normally takes less than 100 ms. For single¬ 
shot EPI, image artifacts due to macroscopic motions are completely eliminated. The advent of 
single-shot, diffusion-weighted EPI has made it possible to obtain quantitative diffusion data in 
vivo from virtually artifact-free images. [80] 

Normally, diffusion-weighted EPI is performed with a spin-echo preparation sequence, and 
a single-shot EPI acquisition. One disadvantage of this approach is that the signal-to-noise in EPI 
images is generally much less than what can be obtained using a conventional, PGSE sequence, 
degrading the precision with which the diffusion measurements can be obtained. The addition of 
difiusion gradients further decreases the signal-to-noise in the images. One possible solution to this 
problem is to reduce the signal loss due to T 2 relaxation by increasing diffusion gradient 
amplitudes, and decreasing the echo time, TE. This would result in a substantial increase in initial 
signal-to-noise at b = 0 s/mm^, aillowing high quality images to be obtained for much higher b 
factors. This approach zu’gues for the use of specialized, high-strength gradient coils to provide 
these higher gradient amplitudes. As an example, we can consider the benefit of increasing the 
gradient strength for each coil to 100 mT/m. Applying these gradients along all three axes gives a 
total gradient strength of 173 mT/m. Using diffusion-pulse timings of 6 = 5 ms, A = 15 ms, the 
smallest TE that can be achieved (using the same acquisition parameters as in our study) is TE = 30 
ms. For Tj = 68 ms (See Appendix 2), the gain in signal achieved at this TE is approximately 50% 
over that obtained in our study with TE = 56 ms. Moreover, a maximum diffusion weighting of b = 
717 s/mm^ could be obtained using these parameters. A secondary benefit to this approach is the 
ability to dramatically decrease the T 2 -weighting in an image. In PGSE images with significant T 2 - 
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weighting, the image contrast is determined by a mix of ADC- and Tj-weighting. Since the ADC 
and Tj reflect different aspects of water mobility within the tissue, these can often have opposing 
effects on signal intensity within a voxel. Further, if several different Tj compartments are present 
within a single voxel, Tj-weighting of these compartments will affect the accuracy of the mean 
ADC measurement for the voxel. 

EPI would enable rapid acquisition of a series of diffusion-weighted images, making it 
technically possible to obtain an entire diffusion decay curve in less than one minute. This would 
provide more information about the character of the diffusion decay than a three-point 
measurement. In addition, the availability of higher-strength gradients would provide stronger 
diffusion-weighting, and increased precision of the diffusion asymptote measurement. According 
to prior work by Turner and Le Bihan, [81] the maximum diffosion-weighting in our study (450 
s/mm^) may not be adequate for complete suppression of the microvascular flow component of the 
signal. In our prior example, where the maximum available gradient strength is increased to 173 
mT/m, the diffusion b factor that could be obtained with TE = 30 ms increased dramatically to b = 
717 s/mml The availability of higher-strength gradients therefore makes it possible to sample the 
asymptotic portion of the decay curve, and thereby determine the molecular diffusion coefficient in 
tissue independent from the microvascular flow. 

5: Development of Hardware for Microvessel/Diffusion Breast Imaging 

A high-strength 3-axis local gradient coil set was constructed for magnetic resonance 
imaging of the breast. This coil design will allow diffusion imaging and high-velocity-resolution 
blood flow imaging around human breast tumours in clinically acceptable exam times. These are 
the first specialized gradient coils for breast imaging to ever be presented (2nd Meeting of the 
Society of Magnetic Resonance, San Francisco, CA, August 1994), and are the first ultra-high- 
strength gradient coils (>200mT/m) designed for human imaging applications (for body parts larger 
than fingers). The coils were tested using a grid phantom, and a diffusion phantom, and 
demonstrated excellent gradient uniformity and accurate diffusion measurement in the region of 
interest that we defined for breast imaging (based on common bra sizes). 

The high gradient efficiencies of these coils will allow much higher spatial-resolution 
imaging than can currently be achieved in a clinical setting, and additionally, will provide very 
large amplitude motion-sensitizing gradients, thus facilitating diffusion imaging or flow imaging in 
small blood vessels while maintaining short echo times. T2 has been measured for human breast 
tissue, with reported values ranging between 60-150 ms. [82,83] Recently, we reported T2 ~ 70 ms 
for human breast tumour (MCF7 human breast carcinoma) growing in a mouse mammary fat pad. 
(Appendix 2) More fibrous tumours can have significantly shorter T2s. For this reason, initial 
signal loss at b=0 s/mm^ due to T2-weighting of the image can severely limit the range of b-factors 
for which images can be obtained, which has the result of limiting the precision of the ADC 
measurement. Our approach is to shorten TE while maintaining diffusion weighting by increasing 
diffusion gradient amplitudes. This gives a substantial increase in initial signal intensity for b=0 
s/mm^, allowing high quality images to be obtained for much higher b-factors. Even in cases where 
signal loss from T2-weighting is not so severe as to limit ADC precision, the ability to obtain high 
diffusion weighting with short TE reduces dramatically the T2 contrast in an image. This will 
allow acquisition of a diffusion-weighted image where the contrast is determined mainly by tissue 
ADC, rather than by a mix of ADC and T2. 

Suggestions for Improvements: 

Imaging results with these coils demonstrated excellent gradient uniformity for two of three 
axes in a large region near the top of the coils, and acceptable gradient uniformity for the third (Y) 
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axis in a smaller region. Initially, the feasibility for using this coil design to measure diffusion and 
microvessel flow in breast lesions will be evaluated clinically by selecting women who have breast 
lesions that can be positioned inside our defined volume of interest; however, improvements on the 
electrical design of the Y gradient coil to provide more uniform gradients over a larger volume 
would allow scanning with a wider range of breast sizes. 

The coils are not torque-balanced, and the horizontal return-wire plates experience strong 
forces in the upward/downward direction. It is possible to constrain the motion of the coils by an 
appropriately designed mounting system, but the design of such a system for human imaging is 
beyond the scope of this thesis, and will not be undertaken immediately. It was therefore decided 
that no human images would be obtained with the initial prototype coils built for this paper, and 
that subsequent prototypes would require extensive testing before an attempt to obtain images of a 
human breast would be made. 

This Paper fincluded as Appendix 3) has been published as : C.F. Maier, H.N. Nikolov, 
K.C. Chu, B.A. Chronik, B.K. Rutt, "Practical Design of a High-Strength Breast Gradient", 
Magnetic Resonance in Medicine, 39(3), pp.392-398. Portions of this work have been presented at: 
the 2nd Meeting of the Society of Magnetic Resonance, San Francisco, CA, August 1994 and the 
3rd Meeting of the Society of Magnetic Resonance, Nice, France, August 1995. This work 
received the University of Western Ontario’s Dean’s A,ward for Excellence in Research (1995), and 
was an invited presentation at the Society of Magnetic Resonance Breast Imaging Workshop, 
Washington D.C., June 24-25,1995. 
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^ Conclusions: 

1. The fibroadenomas in our histological study of lesion vascularity contained, on average, 
more microvessels than some of the invasive carcinomas. This finding directly contradicts the 
assertion by imaging researchers that benign breast lesions are not vascularized. 

2. The numbers of small vessels were higher around the boundaries of invasive carcinomas 
compared to the lesion centers, whereas small vessels were more uniformly distributed in 
fibroadenomas. This characteristic spatial distribution of blood vessels may allow differentiation 
between fibroadenomas and carcinomas using non-invasive imaging methods sensitized to flow 
in microvessels. 

3. Quantitative diffusion measurements were performed in MCF7 human breast tumours 
implanted in nude mice. Apparent Diffusion Coefficients were specific for various components 
of solid tumours (viable tumour, necrosis, fibrous tissue). Diffusion maps were compared to 
histology, and correlated well with gross tumour morphology. These measurements provide a 
foundation for the IVIM experiments that will be undertaken. 

4. PGSE measurements made in a well-vascularized tumour model were highly variable due 
to macroscopic mouse motions during the course of the MR imaging protocol. The precision of 
•these measurements was insufficient for reliable calculation of vascular parameters using IVIM 
theory. These results indicate that higher gradient strengths are necessary to perform fast, motion- 
insensitive diffusion-weighted MR protocols. 

5. An ultra-high-strength 3-axis local gradient coil set was constructed for magnetic 
resonance imaging of the breast. This coil design will allow diffusion imaging and high- 
velocity-resolution blood flow imaging around human breast tumours in clinically acceptable 
exam times. Initial imaging and diffusion measurements using phantoms are promising and 
prove our design concept; however, construction of further prototypes and more extensive testing 
is required before hmnan scanning may be safely undertaken. 
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Figure 1. PGSE pulse sequence for diffixsion/IVIM measurements with TE=75ms, 6=10ms, and 
A=50ms. b-factors for this TE range between 0 and 485 s/mm2. 
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Figure 2. Mean signal intensity inside an ROI defined within a water standard. The acquisition 
time for each image was 4.5 min. 
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Figure 3. Mean signal intensity inside an ROI defined for a mouse tumour. The acquisition time 
for each image was 4.5 min. 
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Invasive Carcinomas and 
Fibroadenomas of the Breast: 
Comparison of Microvessel 
Distributions—Implications 
for Imaging Modalities^ 
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At present, the primary method of the detection and diagnosis of breast disease is 
conventional mammography, which has a proved sensitivity of 85%-95% (1,2). However, 
specificity with, this modality is substantially lower, and additional diagnostic procedures 
are often necessary (3,4). A definitive diagnosis is sometimes possible with noninvasive 
techniques, but, for most lesions, biopsy is necessary to obtain reliable results (5-8). In the 
majority of cases, breast lesions progressing to biopsy are benign (3,9,10). Surgical biopsy m 
these cases has served only to provide a diagnosis, since it is not necessary to remove a 
benign lesion unless it is causing pain (11,12). This has led to the evaluation of iinaging 
methods such as magnetic resonance (MR) imaging and color Doppler ultrasound (US) as 
adjunctive noninvasive diagnostic examinations (13,14). 

Both these modalities rely on the principle that an inaease in the number of blood 
vessels, increased vascular permeability, or both around malignant lesions result from 
tumor angiogenesis and create characteristic, identifiable patterns that can be dishn- 
guished from those associated with benign lesions. Although early contrast material- 
enhanced MR imaging and color Doppler US studies were reported to have extremely high 
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TABLE 1 

Invasive Ductal Breast Carcinomas: Clinical Data and Vessel Counts 


Case 

No./ 

Patient 

Age(y) 

Tumor 

Diameter 

(cm) 

Tumor 

Grade 

Lymph 

Node 

Status 

Vessel Count per Field 

Mean ± Standard Deviation Range 

Ratio of Mean 
in Periphery 
to Mean 
in Center 

. i 

i 

■ 1 

Periphery 

Center 

Periphery 

Center 

1/44 

1.5 

11 

Negative 

13.0 ±2.20 

8.00 ± 2.20 

9-16 

5-12 

i.el ‘ 

■- ! 

2/66 

0.7 

II 

Negative 

20.8 ± 4.40 

11.4 ± 2.10 

16-26 

9-14 

1.82 

■■■ 1 

3/74 

2.0 

1 

Not known 

25.1 ± 6.80 

17.3 ± 3,70 

17-38 

13-23 

1.45 

: 1 

4/47 

1.7 

II 

Positive 

32.0 ± 10.5 

15.6 ± 5.50 

20-53 

8-25 

2.05 

• - • 4 

5/75 

3.0 

11 

Not known 

19.8 ± 3.20 

10.8 ± 2.40 

16-24 

7-15 

1.83 

i 

. z 

6/70 

2.5 

Hi 

Positive 

19.8 ± 2.40 

7.70 ± 1.80 

16-22 

5-10 

2.57 

1 

7/54 

2.5 

il 

Positive 

24.5 ± 4.40 

10.8 ±4.70 

17-30 

4-17 

2.27 

1 

8/55 

1.5 

11 

Negative 

29.1 ± 10.3 

28.8 ± 4.30 

14-49 

23-35 

1.01 


9/49 

2.2 

111 

Positive 

29.0 ± 3.30 

21.1 ± 7.00 

25-32 

14-32 

1.37 


10/30 

1.8 

111 

Positive 

71.0 ± 18.1 

33.9 ± 10.8 

54-105 

24-57 

2.09 


11/80 

2.5 

HI 

Not known 

30.2 ± 6.60 

16,6 ± 3.20 

22^3 

13-23 

1,82 


12/31 

2.0 

III 

Positive 

23.4 ± 4.60 

18,0 ± 4.90 

16-30 

10-27 

1.30 


13/52 

1.3 

11 

Negative 

35.4 ± 13.0 

17.2 ± 3.90 

17-55 

9-22 

2.06 


14/57 

1,0 

1 

Negative 

21.6 ±5.60 

22.8 ± 4.30 

16-31 

18-29 

0.95 


15/79 

2,0 

1 

Not known 

41.6 ±7.50 

35.6 ± 9.50 

28-54 

24-51 

1,17 


16/46 

5.0 

II 

Positive 

23.6 ± 3.00 

19.7 ±2.60 

17-27 

16-24 

1.20 

-.-■r ■ 

17/49 

2.5 

II 

Negative 

30.1 ± 7,60 

18.2 ± 3.30 

20-45 

13-22 

1.65 


18/47 

1.5 

1 

Negative 

52.5 ± 14.9 

42.5 ± 6.30 

39-88 

36-57 

1.24 


19/61 

0.4 

11 

Not known 

55.1 ± 12.8 

25.9 ± 6.40 

35-72 

18^0 

2.13 



diagnostic specificity (15-17), further 
studies with larger sample sizes did not 
achieve the same accurate delineation of 
breast carcinomas. In particular, a substan¬ 
tial fraction of fibroadenomas were indis¬ 
tinguishable from malignant lesions (13,18- 
20). The. inability to accurately identify 
fibroadenomas is a considerable flaw for a 
potential adjunctive modality, because 
these lesions often resemble solid carcino¬ 
mas at both mammography and clinical 
presentation (21,22) and account for 30%- 
50% of biopsies of benign lesions (23). 

Recently, several MR groups have ob¬ 
served a pattern of heightened contrast 
material enhancement around the periph¬ 
ery of a lesion in 15%-30% of malignant 
lesions (24,25). This ''positive rim sign" 
(18) has also been observed by using color 
Doppler US. It has been suggested that 
this criterion could aid in separating ma¬ 
lignant from benign lesions. 

Several researchers (25-27) have also 
compared MR contrast material enhance¬ 
ment parameters with numbers of blood 
vessels observed at histologic analysis. In 
these studies, vessels were counted in 
areas of highest vascularity or "hot spots," 
potentially valuable ^prognostic indica¬ 
tors for invasive carcinomas (28). Two of 
the three researchers (25,26) found a sta¬ 
tistically significant correlation between 
the degree of enhancement and microves¬ 
sel density in hot spots. Unfortunately, 
researchers in none of these studies at¬ 
tempted to evaluate the spatial distribu¬ 
tion of vascularity over entire lesions. In 
addition, no data are available on the 


comparison of color Doppler results with 
histologic confirmation of the presence 
of larger vessels. 

To use vascular imaging techniques ef¬ 
fectively, it is necessary to first identify a 
difference in the vascular pattern be¬ 
tween malignant and benign lesions that 
is prominent enough to allow differentia¬ 
tion of these lesions. In addition, the 
study of the differences, similarities, or 
both between vascular distribution and 
density in malignant and benign breast 
lesions could serve to illuminate the 
unique mechanisms contributing to the 
formation of a microcirculatory network 
in a malignant tumor and its subsequent 
role in tumor metastasis. With these objec¬ 
tives in mind, we have evaluated spatial 
patterns of blood vessel distribution and 
quantitatively assessed microvessel den¬ 
sity in invasive ductal carcinomas and 
fibroadenomas. 


MATERIALS AND METHODS 

Routinely processed, formalin-fixed, par¬ 
affin-embedded tissue from'39 lumpecto¬ 
mies in 19 patients with 19 invasive 
ductal carcinomas and in 20 patients 
with 20 fibroadenomas was retrieved from 
the files of the Department of Pathology 
at London Health Sciences Centre, Lon¬ 
don, Ontario, Canada. Patients with carci¬ 
noma ranged in age from 30 to 80 years 
and patients with fibroadenomas ranged 
in age from 19 to 49 years. Immunohisto- 
chemical studies to detect factor VIII- 


related antigen (rabbit polyclonal anti- 
body; Dimension Labs, Mississauga, Ontario, 
Canada; 1:100 dilution) were performed 
on 4 jjLm-thick recut sections by using the 
avidin-biotin peroxidase technique. All > 
tissues were counterstained with hema¬ 
toxylin. A standard hematoxylin-eosin 
slide was also prepared from an adjacent 
slice of the fibroadenomas for classifica¬ 
tion purposes. 

Fibroadenomas were classified by one 
of the investigators (M.M.) by using the . 
three categories defined by Ohtani and 
Sasano (29)—myxoid, fibrocellular, and 
sclerosed—^with the addition of two cat¬ 
egories: mixed lesions containing both 
myxoid and fibrous characteristics (fibro- 
myxoid) and lesions showing only cellu¬ 
lar characteristics (cellular). Patient age, 
tumor diameter, tumor grade (Scarff-Bloom- 
Richardson criteria) (30), and lymph node 
status were obtained from the surgical 
pathology consultation report. 

To evaluate the spatial distribution of 
capillaries and smaller arterioles and ven¬ 
ules, microvessels smaller than 40 \xm in 
diameter (31) were counted by using light 
microscopy on X200 fields (x20 objec¬ 
tive lens, xlO ocular lens) by one of the 
investigators (M.M.). Vessel identifica¬ 
tion was performed according to the crite¬ 
ria used by Weidner et al (28). Any brown- 
staining, freestanding endothelial cell or 
group of endothelial cells was counted as 
a vessel. 

Ten fields were chosen at random in 
both the central and peripheral regions of 
each tumor (provided there was sufficient 
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pibroadenomas; Clinical Data and Vessel Counts 

Case 

No./ 

patient 

Age(y) 

Tumor 

Diameter 

(cm) 

Tumor 

Classification 



Vessel Count per Field 



Ratio of Mean 
in Periphery 
to Mean 
in Center 

Mean 

± Standard Deviation 


Range 


Periphery 

Center 

Norma! 

Periphery 

Center 

Normal 

Tomi 

1.0 

Fibromyxoid 

34.9 ±13.2 

25.3 ± 5.68 

43.6 ± 9.02 

15-57 

17-35 

32-55 

1.42 

71/47 

1.4 

Fibromyxoid 

33.1 ± 7.34 

36.7 ± 10.0 

79.6 + 16.5 

24-44 

19-58 

54-92 

0.9 

77/43 

1.0 

Cellular 

112 + 29.8 

95.6 ± 24.4 

133 ± 58.2 

54-157 

76-137 

75-200 

1.17 

73/19 

3.5 

Cellular 

35.2 ±11.9 

17.0 ±4.81 

55.0 ± 6.82 

20-57 

7-25 

43-60 

2.07 

74/44 

1.0 

Fibrocellular 

48.5 ± 13.0 

53.4 ± 10.4 

75.6 ± 13.3 

29-68 

40-67 

60-96 

0.91 

75/^7 

0.6 

Sclerosed 

47,9 ± 9.75 

40.8 ± 5.26 

46.8 ± 5.93 

31-65 

32-46 

37-52 

1.17 

7A/38 

0.6 

Myxoid 

51.7 ± 9.04 

68.0 ± 19.8 

98.2 ± 32.1 

41-70 

54-82 

65-151 

0.76 

77/43 

2.0 

Cellular 

39.3 ± 13.1 

45.7 ± 12.6 

no ± 26.2 

20-60 

25-71 

85-144 

0.86 

78/31 

1.5 

Sclerosed 

39.1 ± 7,23 

27.3 ± 8.96 

75.7 ± 16,6 

26-47 

14-39 

57-89 

1.43 

29/39 

2.5 

Fibrocellular 

29.1 ±11.0 

18.0 ±6.02 

39.6 ± 8.85 

15-50 

7-27 

28-51 

1.62 

30/49 

1.0 

Myxoid 

40.4 ±11.0 

24.0 ± 9.90 

99.2 ± 58.5 

27-58 

17-31 

39-167 

1.68 

31/22 

2,7 

Cellular 

24.7 ± 6.93 

15.0 ±0.0 

53.2 ± 37.2 

14-34 

15 

21-113 

1.64 

32/34 

1.7 

Fibromyxoid 

30.7 ± 15.8 

34.4 ± 12.3 

59.2 ± 12.0 

14-59 

11-51 

45-71 

• 0.89 

33/35 

2.0 

Sclerosed 

37.4 ± 8.28 

30.6 ± 14.1 

0 

22-50 

15-58 

0 

1.22 

34/21 

1.1 

Fibromyxoid 

64.9 ±12.0 

64.6 ±12.2 

^ 0 

47-78 

46-88 

0 

hO 

35/27 

1.4 

Myxoid 

51.2 ± 8.89 

45.8 ± 10.4 

100 ±49.2 

41-65 

29-63 

54-152 

1,12 

36/37 

1.2 

Myxoid 

50.8 ± 14.9 

52.6 ± 11.5 

50.5 ± 24.8 

33-77 

34-71 

33-68 

0.97 

37/20 

2.5 

Fibromyxoid 

37.0 ± 9.57 

21.8 ± 8,47 

96.8 ± 25.0 

23-55 

11^0 

66-130 

1.7 

38/37 

1.2 

Myxoid 

46.0 ± 18.0 

30.8 ± 6.58 

86.8 ± 19.8 

23-75 

19-38 

55-102 

1.49 

39/47 

1.2 

Myxoid 

65.0 ± 14.4 

66.2 ± 13.9 

146 ±77.9 

48-90 

48-81 

66-229 

0.98 


tissue). A peripheral field was defined as a 
field inside the tumor boundary (ie, an 
imaginary line connecting the most pe¬ 
ripheral nests was considered the leading 
edge of invading cells in the case of invasive 
carcinomas) with the microscopic field 
circumference touching the boundary, A 
central field was defined as a microscopic 
field inside the tumor at least one x200 
field of view away from the boundary. 
Central fields were not defined in ne¬ 
crotic, avascular regions of the carcino¬ 
mas; however, tumors containing necro¬ 
sis were noted. For fibroadenomas, up to 
five fields of normal tissue adjacent to the 
lesion were also evaluated (provided there 
was sufficient normal tissue within the 
slide limits): 

No tissue present on a carcinoma slide, 
regardless of histologic appearance, was 
considered normal (normal-appearing 
lobules from cancerous breasts have 
shown an angiogenic up-regulation when 
compared with lobules from noncancer- 
ous breasts [32]). Microvessels 40 \xm or 
larger in diameter were also qualitatively 
evaluated on each slide. 

Student t tests were performed for indi¬ 
vidual lesions to determine whether there 
was a statistically significant difference 
between the vessel density in the periph¬ 
ery and that in the lesion center. (A 
Student t test could not be performed on 
one fibroadenoma [case 31]; size permit¬ 
ted only one countable central field ac¬ 
cording to our criteria.) Ratios of mean 
peripheral to mean central vessel counts 


were calculated for each lesion, and t tests 
were performed to assess differences in 
vessel coxmts or in this ratio between inva¬ 
sive carcinomas and fibroadenomas and 
between carcinomas with lymph node 
involvement and carcinomas without 
lymph node involvement. One-way analy¬ 
ses of variance were performed to evalu¬ 
ate differences in vessel counts between 
the five categories of fibroadenomas and 
also between the three grades of invasive 
carcinomas (Scarff-Bloom-Richardson cri¬ 
teria [30]). In all cases, a significant differ¬ 
ence was defined at P less than ,01. 


RESULTS 


Mean vessel counts (vessel density within 
a X200 field) and their standard devia¬ 
tions and ranges are shown in Tables 1 
and 2 for invasive ductal carcinomas and 
fibroadenomas, respectively. Invasive 
ductal carcinomas ranged in diameter 
from 0.4 to 5.0 cm, while fibroadenomas 
ranged in diameter from 0.6 to 3.5 cm. The 
invasive ductal carcinomas could be cat¬ 
egorized as follows: four (21%) grade I 
carcinomas, two without lymph node 
involvement and two in which no lymph 
node dissection was performed; 10 (53%) 
grade II carcinomas, five without lymph 
node involvement, three with lymph 
node involvement, and two in which no 
lymph node dissection was performed; 
and five (26%) grade III carcinomas, all 
with lymph node involvement. Classi¬ 


fication of the fibroadenomas yielded six 
(30%) myxoid, five (25%) fibromyxoid, 
four (20%) cellular, two (10%) fibrocellu- 
lar, and three (15%) sclerosed lesions. 
Large necrotic centers (ie, >30%-40% of 
the lesion was necrotic) were evident in 
two (10%) of 19 carcinomas (cases 7 and 
10). In one of these tumors (case 10), the 
area of necrosis was surrounded by a 
prominent ring of small blood vessels. 
Smaller areas of partial necrosis (approxi¬ 
mately 10%~15% of the lesion was ne¬ 
crotic) were observed within the lesion 
center in an additional eight (42%) of the 

19 invasive ductal carcinomas. 

The number of vessels in a X200 field 
for the invasive ductal carcinomas ranged 
between 9 and 105 (mean, 31.4) for the 
peripheral and 4 and 57 (mean, 20.1) for 
the central fields. In previous vessel den¬ 
sity studies, only areas of highest vascular¬ 
ization within the tumor (tumor hot 
spots) were selected for vessel counting 
(28,33,34). Maximum vessel counts for 
invasive ductal carcinomas in our study 
were within these reported ranges. For 
fibroadenomas, the range of vessel densi¬ 
ties for the peripheral fields was 14-157 
vessels (mean, 46,0 vessels) and that for 
the central fields was 7-137 vessels (mean, 
40.7 vessels). Normal tissue surrounding 
the fibroadenomas (present in 18 of the 

20 cases) had a range of 21-229 vessels 
(mean, 80.5 vessels) per field. In 16 (89%) 
of the 18 cases, vessel counts in normal 
tissue were significantly higher than those 
in fibroadenomas. It is interesting to note 
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that the highest normal tissue vessel 
counts were found in three women aged 
43^9 years, while the mean age of pa¬ 
tients with fibroadenomas was 35 years. 
For these three women, the normal breast 
tissue was noticeably more vascular than 
that in other patients and was predomi¬ 
nantly composed of dense lobular and 
ductal tissue, with little fatty tissue pres¬ 
ent. 

Vessel counts for invasive ductal carci¬ 
nomas are shown in Figure la; those for 
fibroadenomas are shown in Figure lb 
(cases have been arranged on the x axis in 
order of increasing tumor diameter). Mi¬ 
crovessel density was significantly higher 
(P < .01) in the tumor periphery than it 
was in the tumor center in 15 (79%) of 
the 19 invasive ductal carcinomas and in 



ca ca ca ca ca ca ca ca ca ca ca ca ca ca ca ca ca ca ca 

19 2 14 13 1 18 8 4 10 3 12 15 9 17 7 6 11 5 16 


Carcinoma Cases 




three (16%) of 19 fibroadenomas. 

Figure 2 shows high magnification 
fields in the periphery (Fig 2a) and center 
(Fig 2b) of a typical invasive ductal carci¬ 
noma (case 13). In one of the four carcino¬ 
mas in which no significant difference 
was found between peripheral and cen¬ 
tral vessel densities (case 18), the carci¬ 
noma consisted Of many nests of carci¬ 
noma cells with distinctive rings of small 
vessels surrounding each of the nests. The 
mean ratio of peripheral to central vessel 
counts was 1.66 ± 0.46 (standard devia¬ 
tion; range, 0.95-2.57) and 1.25 ± 0.36 
(range, 0.76-2.07) in invasive ductal carci¬ 
nomas and fibroadenomas, respectively, 
and the ratios of peripheral to central 
vessel counts was significantly different 
between these groups (P < .01). 

The individual peripheral versus cen¬ 
tral vessel density ratios for both groups 
of lesions are compared in Figure 3. No 
correlation was observed between mean 
peripheral counts, mean central counts, 
maximum peripheral counts, maximum 
central counts, or in the peripheral versus 



Fibroadenoma Cases 


h- 

Figure 1. Mean vessel counts of lesions arranged in order of increasing tumor diameter. * - those 
cases in which the vessel density in the periphery was significantly greater than that in the center 
(? < .01). (a) Graph shows mean vessel counts in the peripheral and central fields of invasive 
ductal carcinomas {ca), (b) Graph shows mean vessel counts in the peripheral and central fields ot 
fibroadenomas (fa) and mean vessel counts in the surrounding normal tissue. + - the fibroadenoma 
case in which only one central field was present and no statistical analysis could be j 

Fibroadenoma cases 1-20 correlate with cases 20-39 in Table 2 (eg, fa I - case 20 in Table 2, fa 2 
case 21). 


central vessel density ratios with carci¬ 
noma grade or lymph node status. 

One-way analyses of variance per¬ 
formed between the maximum periph¬ 
eral, maximum central, mean peripheral, 
and mean central counts of the fibroad¬ 
enomas yielded no statistically signifi¬ 
cant difference (P < .01) in vessel densi¬ 
ties between the five classifications. All 
three sclerosed fibroadenomas in our 
study had numerous blood vessels, com¬ 
parable even to the number of vessels 
observed for the cellular and myxoid 
groups. An example of a sclerosed but 
highly vascularized fibroadenoma (case 
28) is shown in Figure 4, and a typical 
cellular fibroadenoma (case 22) is shown 
in Figure 5 for comparison. There was 
also no correlation of any vascular param¬ 


eters with patient age or tumor size for 
either fibroadenomas or carcinomas. 

In eight (42%) of the 19 invasive ductal 
carcinomas, vessels 40 |xm or larger in 
diameter were observed just outside the 
leading edge of the tumor (Fig 6). In four 
(21%) of these cases, larger vessels could 
also be seen immediately inside the invad¬ 
ing boundary. In five additional carcino¬ 
mas (26%), a patch of larger vessels was 
observed within a peripheral field of a 
tumor. Larger vessels were not observed 
within the centers of any of the carcino¬ 
mas in our study. In 10 (50%) of the 20 
fibroadenomas, the presence of larger ves¬ 
sels was also noted. In six of these cases 
(30%), there were several large vessels 
present in the periphery of the lesion; 


while in two cases (10%), larger vessels in 
the periphery were abundant. In the re¬ 
maining two fibroadenomas (10%), many 
large vessels were uniformly distributed 
throughout the lesion. These vessels were 
easily observed with low magnification. 


DISCUSSION __ 

We identified statistically significant over¬ 
lap in microvessel counts between the 
invasive ductal carcinomas and fibroad¬ 
enomas in our study. This finding may 
help explain the overlap in enhancement 
characteristics between invasive breast 
carcinomas and some fibroadenomas on 
contrast-enhanced MR imaging studies. 
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Figure 2. Case 13. T 5 ^ical invasive ductal carcinoma, (a, b) Pathologic slices show carcinoma in 
the (a) periphery and (b) center. (Stained for factor Vlll-related antigen and counterstained with 
hematoxylin; original magnification, x200.) In a and b, arrows show vessels outlined in brown 
(factor Vlll-related antigen causes endothelial cells that line the lumen of the blood vessels to 
stain brown). 
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Figure 3. Histogram shows the distribution of mean peripheral versus mean central vessel count 
ratios for invasive ductal carcinomas and fibroadenomas. 


Enhancement characteristics on a con¬ 
trast-enhanced MR image reflect blood 
volume (microvessel density), blood flow, 
and vessel permeability, which are all 
vasculature characteristics that have been 
linked to tumor angiogenesis. The ob¬ 
served overlap in MR imaging enhance¬ 
ment curves and in microvessel density 
between carcinomas and fibroadenomas 
implies that the same processes may be 
controlling the vascular development in 
both of these classes of lesions. 

It is possible that the numerous small 


vessels we observed in fibroadenomas 
may be the result of angiogenic stimula¬ 
tion initiated by the fibroadenoma. This 
hypothesis is also supported by the in¬ 
creased number of vessels 40 jxm or larger 
observed in the periphery of the fibroad¬ 
enomas, similar to the pattern of periph¬ 
eral blood vessels known to be associated 
with angiogenesis in carcinomas (35). It is 
well known that premalignant lesions 
can induce angiogenesis, and it has been 
shown that benign breast tissue removed 
from a location distal to a carcinoma but 


within the same breast can induce new 
blood vessel formation in a rabbit comeal 
model (32,36). Therefore, it is not unreal¬ 
istic to suggest that other benign breast 
diseases may also achieve an angiogenic 
state. 

Our results support the observation 
that microvessel density is higher in the 
active tumor periphery of invasive ductal 
carcinomas than in the tumor center. 
These findings are consistent with con¬ 
trast-enhanced MR imaging findings of 
rim enhancement in carcinomas. In the 
study by Buadu et al (25), rim-enhancing 
carcinomas were observed to have vary¬ 
ing degrees of central desmoplasia, associ-' 
ated with lower vessel density. In the 
work of Orel et al (24), rim enhancement 
was observed in five of 16 carcinomas, 
but none of the rim-enhancing carcino¬ 
mas exhibited central necrosis. 

In our group of lesions, a significantly 
higher microvessel density was observed 
in the lesion periphery in 15 (79%) of the 
19 invasive ductal carcinomas. Only two 
of these cases (cases 7 and 10) had well- 
defined necrotic cores, and, moreover, 
central fields were not defined within 
these avascular regions; rather, they were 
limited to areas of viable tumor. An addi¬ 
tional eight cases had small necrotic re¬ 
gions within the tumor; however, only 
five of these had significantly higher pe¬ 
ripheral vessel counts than central vessel 
counts. Although the mean microvessel 
density was significantly higher in periph¬ 
eral fields than in central fields in most 
carcinomas in our sample, there was still 
considerable overlap between peripheral 
field counts and central field counts 
within each lesion. This may explain why 
the percentage of carcinomas reported to 
exhibit rim enhancement on MR images 
is much smaller than the percentage of 
lesions we found to have increased periph¬ 
eral vessel density (24,25). 

Vessel density in fibroadenomas was 
more uniform throughout the tumor than 
it was in carcinomas, with no statistically 
significant difference between periphery 
and center in 16 (84%) of 19 cases. Areas 
of increased vessel density did exist within 
fibroadenomas; however, these areas usu¬ 
ally corresponded to a single or several 
lobules, within which a uniformly ele¬ 
vated distribution of blood vessels was 
observed. These results are in accord with 
the homogeneous MR imaging contrast- 
enhancement patterns observed for fibro¬ 
adenomas and are suggestive of possible 
improvements in specificity that may re¬ 
sult from the addition of the evaluation 
of morphologic features to an MR imag¬ 
ing examination (25,37,38). 
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Two authors (24,39) have compared con¬ 
trast-enhancement characteristics with his¬ 
tologic findings in an attempt to explain the 
rapid, marked contrast material enhance¬ 
ment observed in some fibroadenomas. 
In these studies (24,39), nonsclerosed fibro¬ 
adenomas enhanced strongly, but the en¬ 
hancement of sclerosed fibroadenomas 
was highly variable, with some sclerotic 
lesions also exhibiting statistically signifi¬ 
cant rapid signal intensity increases. This 
finding was contrary to the expectation 
that sclerotic lesions, being acellular, 
would have fewer blood vessels than myx¬ 
oid or cellular lesions. In the group of 
fibroadenomas in our study, no statisti¬ 
cally significant decrease in the number 
of microvessels was observed in the three 
sclerotic lesions as compared to the num¬ 
ber of microvessels in the cellular or 
myxoid category. Lesions in all these 
groups contained areas of low microves¬ 
sel density and areas of high microvessel 
density, with no general trend observable 
with category; however, the number of 
fibroadenomas in each of our categories 
was small, and further investigation is 



a. 

Figure 4. Case 28. Pathologic sUces show a sclerosed but highly vascularized fibroadenoma^ 
stained for factor Vlll-related antigen and counterstained with hematoxylin, (a) Pathologic slice^ 
at X40 magnification, (b) Part of a X200 field in a pathologic slice. Salient features of this| 
fibroadenoma stage are few stromal cells (c in b), hyalinized collagen fibers (h in b), and^ 
compressed ducts (d in a). > j 


required before any definite conclusions 
can be drawn. 

Our observation of larger microvessels 
in the periphery of both carcinomas and 
fibroadenomas is consistent with the color 
Doppler US study by Kubek et al (18), in 
which a preponderance of detected blood 
flow near the peripheries of both malig¬ 
nant lesions and fibroadenomas was re¬ 
ported. This finding is not inconsistent 
with the fact that MR imaging rim en¬ 
hancement is observed only in malignant 
lesions, since color Doppler US is sensi¬ 
tive to much larger vessels than is con¬ 
trast-enhanced MR imaging, and our find¬ 
ings show that the spatial distribution of 
these larger vessels is different from that 
of smaller vessels. Two investigators 
(18,19) have reported that the number of 
vessels observed by using color Doppler 
US in fibroadenomas is increased with 
increasing lesion size. We have not noted 
any correlation between the number of 
vessels 40 iJim or larger and fibroadenoma 
size. Very few vessels larger than 200 |xm 
could be observed on a single slide; it was 
therefore not possible to conclude any¬ 
thing about the' spatial distributions of 
Doppler-detectable vessels in our study. 

Numerous investigators have estab¬ 
lished microvessel density in hot spots as 
an important prognostic indicator in 
breast carcinoma (28,40). Several authors 
have compared MR contrast-enhance¬ 
ment patterns with vessel density in hot 
spots; these investigations have been 
more directly relevant to assessing the 



Figure 5. Case 22. Pathologic slices show a typical cellular fibroadenoma stained for facto 
Vlll-related antigen and counterstained with hematoxylin, (a) Pathologic slice at x40 magmnca 
tion. (b) Part of a X200 field in a pathologic slice. Salient features of this fibroadenoma stage ar , 
high numbers of stromal cells (c in b) and open ducts {d in a). 


applicability of MR to predicting disease tions. While most recent work in which 

outcome rather than to determining US has been used to image breast carci- 

whether a lesion is malignant or benign, noma has focused on color Doppler 

We have evaluated microvessel density substantial improvements in vessel detecP 

over entire tumors to provide more perti- ability are expected with the use of intra- 

nent information for using these noninva- venous US contrast material, and it V 

sive techniques for diagnostic applica- possible that sensitivity will be increase 
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Figure 6. Case 7. Pathologic slice of an inva¬ 
sive carcinoma at X40 magnification (stained 
for factor Vlll-related antigen and counter- 
stained with hematoxylin). Note the abun¬ 
dance of larger vessels (>200 ij-m; arrows) 
around the periphery. 

to a level at which capillaries and other 
small microvessels will be visualized. Fur¬ 
ther gains in spatial resolution and vessel 
detectability with Doppler US and MR 
imaging will increase the potential for 
using these imaging modalities to iden¬ 
tify tumor hot spots and noninvasively 
measure microvessel density. 

While further histologic studies are nec¬ 
essary to evaluate vascular patterns in 
nonpalpable lesions, such as in atypical 
hyperplasia and carcinoma in situ, the 
results reported here suggest that the 
additional information on the spatial dis¬ 
tribution of vessels and the ability to 
resolve vessel size may add valuable infor¬ 
mation to the diagnostic work-up of an 
indeterminate solid lesion. These exami¬ 
nations have important implications for 
younger women with palpable lesions, 
since diagnosis by using mammography 
is particularly unreliable in this age group 
and fibroadenomas are common. 

Acknowledgment: The authors thank Elaine 
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Quantitative Diffusion Imaging in Implanted Human 
Breast Tumors 

Cynthia F. Maier, Yael Paran, Peter Bendel, Brian K. Rutt, Hadassa Degani 


Quantitative diffusion measurements were performed in tu¬ 
mors arising from inoculation of nude mice with two human 
breast cancer cell lines (MCF7 and T47D) to evaluate the 
specificity of this technique for characterizing solid tumors. 
ADC maps were compared to histology and correlated well 
with gross tumor morphology. Measured ADCs were highly 
specific for viable and necrotic tumor in the five T47D tumors 
included in this study (P < 0.02), while only two of the five 
MCF7 tumors studied developed distinguishable areas of ne¬ 
crosis. No statistically significant difference was observed in 
ADCs from viable tumor between the different cell lines (P > 
0.05). 

Key words: magnetic resonance imaging; diffusion; breast 
cancer; necrosis. 

INTRODUCTION 

Magnetic resonance diffusion imaging is a technique that 
has recently been applied to imaging human tissues in 
vivo. The rate of water diffusion in tissue directly reflects 
intrinsic tissue parameters such as tissue microstructure 
and the degree of water compartmentalization, and the 
concentration of proteins and macromolecules in the wa¬ 
ter’s environment. In tumors, water diffusion is an im¬ 
portant process that has implications for nutrient deliv¬ 
ery and clearance. 

Preliminary diffusion measurements in animal tumor 
models using MRI have demonstrated potential for pro¬ 
viding information about gross tumor morphology and 
for detecting tumor necrosis (1-3). It has been suggested 
that tumor response to chemotherapy or radiation ther¬ 
apy could be monitored non-invasively using an MR- 
imaging technique based on diffusion (3). To date, how¬ 
ever, most spatially resolved, in vivo studies of water 
diffusion have been based on a single diffusion-weighted 
image, which provides only qualitative information. A 
more quantitative approach is necessary, where apparent 
diffusion coefficients (ADCs) are measured with high 
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spatial resolution for the tumor and the specificity of 
these ADCs for characterizing tumors is assessed. In re¬ 
cent work by Eis et al (4), quantitative analysis of the 
ADC in rat brain tumors demonstrated potential for dif¬ 
ferentiating tumors from normal brain tissue. However, 
the low spatial resolution of images in this study did not 
allow differentiation of morphologically distinct features 
(e.g., areas of necrosis or fibrosis) within the tumors. 

The purpose of this work was to demonstrate the po¬ 
tential for diffusion imaging as a method for characteriz¬ 
ing solid tumors, and for providing information about 
tumor morphology on a cellular level. We evaluated wa¬ 
ter diffusion by using MRI in two different mouse tumor 
models of human breast cancer (MCF7 and T47D cell 
lines). A range of tumor sizes were studied, where some 
of the tumors were allowed to grow to large sizes, result¬ 
ing in cell death and necrosis in the interior of the tumor. 
ADC maps were calculated for the tumors, and average 
ADCs measured for regions of viable neoplastic tissue 
and for necrosis. These quantitative diffusion measure¬ 
ments were compared to histopathology and the speci¬ 
ficity of the ADC for identifying different tumor compo¬ 
nents was evaluated. 

MATERIALS AND METHODS 

Tumor Model 

MCF7 and T47D human breast cancer cells were cultured 
and prepared for inoculation as described previously (5, 
6). Ten female CDl-NU athymic mice (8 weeks old), were 
subcutaneously inoculated in the mammary fat pad with 
approximately 10^ cells in a volume of 0.4 ml (five mice 
for each tumor model). Before injection of the cancer 
cells, a pellet of 17/3 estradiol was implanted under the 
skin of the back (0.72 mg/pellet, 60 day release. Innova¬ 
tive Research of America, Sarasota, FL) to facilitate tumor 
growth. The tumors were allowed to develop for 4-12 
weeks, to sizes ranging between 0,5 and 2.0 cm in diam¬ 
eter. In larger tumors, regions of necrosis developed 
mainly in the centers of the tumors, presumably as a 
result of inadequate nutrient and oxygen delivery to 
these areas. 

MRI Measurements 

Mice were anesthetized before imaging by using Rental 
diluted in PBS solution (Pentobarbitone sodium, CTS 
Chemical Industries, Tel-Aviv, Israel) by an intraperito- 
neal (IP) injection (0.06 mg/g). Follow-up injections (0.03 
mg/g) were administered IP on an hourly basis after the 
initial injection, via an infusion line that extended out¬ 
side the imager. To prevent motion during the examina¬ 
tion, the mice were immobilized for the course of the MR 
experiments by laying them in prone position on a home- 
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ode mouse tray and taping the lower half of the body 
urely to the tray. The upper abdomen and chest were 
ot constrained so that respiration would not be im- 
" jred. A water standard was placed near the tumor for 
ach of the experiments. The room temperature was 
Igintained during the examination at 25-30°C. 

MR images of the tumors were collected at 4.7 Tesla 
bv using a horizontal bore Biospec 4.7/30 animal imager 
(Bruker, Karlsruhe, Germany), and a custom-built 7.5-cm 
diameter radiofrequency (RF) coil. Gradient hardware 
consisted of unshielded gradient coils with a maximum 
gradient strength of 48 mT/m at a rise time of 500 ps, 
using standard gradient pre-emphasis provided by the 
manufacturer. A transverse pilot scan was used to define 
coronal slices covering the whole tumor. A standard 
multislice PGSE sequence was used with TE/TR = 48 
ms/2400 ms, A = 22.4 ms, 8 = 16.2 ms (A and 8 are 
defined as in Fig. 1). TE was kept as short as possible to 
minimize signal loss due to relaxation (our previous 
measurements of in these tumors gave an average 
value of approximately 68 ms). A series of 13 diffusion- 
weighted images was obtained, with diffusion b values 
ranging from b = 0 s/mm^ to b = 754 s/mm^ (b = 
where y is the magnetogyric ratio and g is 
the diffusion-encoding gradient strength as defined in 
Fig. !)• Nine of these 13 images were sensitized to inco¬ 
herent blood flow with b values less than 220 s/mm 
(7-10). Diffusion-weighting was varied by changing the 
amplitude of the diffusion gradients while keeping all 
timing parameters constant. Phase-cycling was per¬ 
formed to cancel the effects of RF inhomopneities on the 
180° refocusing pulse. Images were obtained with 3 cm 
field of view, slice thickness of 1 mm, a matrix size of 
64 X 256 pixels (in-plane resolution of 472 X 118 pm) 
and two averages for an acquisition time of 5 min. and 
10 s. The raw data were zero-filled and Fourier trans¬ 
formed to give images with 256 X 256 pixels. 

For the imaging gradient strengths required at this 
spatial resolution, significant systematic errors are intro¬ 
duced into the determination of diffusion coefficients by 
the interaction of the imaging gradients with the di^- 
sion gradients (11,12). We employed a method for elim¬ 
inating this contribution to the b weighting that was first 
proposed by Neeman et al. (13), where images are ac¬ 
quired by using both positive and negative diffusion 
gradients for each b value and the geometric average of 
these images is calculated. This method also allows the 
elimination of unwanted contributions to the b value 
from the interaction of the diffusion gradients with back¬ 
ground gradients arising from susceptibility variations 
within the sample. 

Following the diffusion imaging sequence, high-reso- 
lution Ta-weighted spin echo images were acquired for 
the same slice locations to provide a reference for the 
diffusion data. Such Ta-weighted images have previously 
been correlated with histology (5) and allow identifica¬ 
tion of necrotic regions, viable tumor, and fibrous tissue. 
The images were acquired using TE/TR = 68 ms/2400 
ms, matrix size = 128 X 256 pixels (in-plane resolution 
of 236 X 118 pm) and four averages. The other parame¬ 
ters were as for the diffusion images. The acquisition 
time for the Tz-weighted images was 20 min. The total 


length of time required for the MR examination including 
setup and localizers was approximately 3 h per animal. 

Image Analysis 

Images were transferred to a DEC 3000 workstation for 
analysis where image analysis was performed by using 
MATLAB software Version 4.2 (The Mathworks, Inc. 
Natick, MA). One image for each b value was obtained by 
geometrically averaging positive-gradient and negative- 
gradient images according to the method proposed by 
Neeman (13). 

Average ADCs were calculated for regions of interest 
(ROIs) in areas that were identified as viable or as ne¬ 
crotic by comparison to the high-resolution Tj-weighted 
scans. These regions were selected on the T 2 -weighted 
scan based on their relative signal intensities, and mor¬ 
phological appearance. The regions were selected as be¬ 
ing “clearly” identifiable as either necrotic or viable. 
Where possible, one ROI for each of viable and necrotic 
tissue were chosen per slice. For these calculations, the 
standard deviation of signal intensities within the ROI 
was taken as the uncertainty in the mean signal intensity. 
Fits having r < 0.97 (14) were rejected. The ADC of the 
water standard was obtained as a reference for each tu¬ 
mor. 

Maps of the apparent diffusion coefficient were calcu¬ 
lated for each slice on a pixel-by-pixel basis from the 
slope of a linear least-mean-squares fit of the natural 
logarithm of signal intensity to the b factor. The noise in 
signal intensity for each pixel was estimated by the stan¬ 
dard deviation of signal intensities in a ROI in the back¬ 
ground. 

Histology 

Immediately after the MR imaging experiments, two mice 
with MCF7 tumors and two mice with T47D tumors were 
sacrificed to provide a comparison to histology. The tu¬ 
mors were surgically removed, fixed in 10% formalde¬ 
hyde solution, imbedded in paraffin, and sectioned to 
correspond to the MR imaging planes. Representative 
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FIG. 1. Multislice PGSE sequence used for diffusion measure¬ 
ments. TE/TR ^or the sequence = 48 ms/2400 ms. A = 22.4 ms, 
8 = 16.2 ms. Diffusion-encoding gradients were on the slice- 
select axis. 
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FIG. 2. Representative MR images from implanted breast tumors, (a) Ta’Weighted image of a 
4.week-old MCF7 tumor (M3) with 128 phase encodes, (b) Diffusion-weighted image of tumor M3 
with b = 636 s/mm^ and 64 phase encodes, (c) A fg-weighted image of a 4-week-o(d T47D tumor 
0*5) with 128 phase encodes, (d) Diffusion-weighted image of tumor T5 with b = 636 s/mm^ and 
64 phase encodes. 


4-/xm thick slices from each section were stained with 
Hematoxylin-Eosin and by using a modified trichrome 
technique to stain mucopolysaccharides of fibrous tissue. 

RESULTS 

In previous work in our laboratory (5,15—17), we studied 
a tumor model where breast cancer cells were implanted 
in the mouse flank. Tumors from both cell lines that were 
studied in this work developed faster in the mammary fat 
pad than had previously been observed in the flank. 
Average ages and diameters for the 10 tumors included in 
this study are given in Table 1. After 4 weeks, most of the 
tumors had reached a size of 0.5-1.0 cm in diameter. In 6- 
to 8-week-old tumors, it was possible to see regions of 
necrosis, usually near the center of the tumor. These 
could be observed both on the T 2 -weighted images, and 
on the highly diffusion-weighted images. After 11 weeks, 
one of the MCF7 tumors had reached a size of 2.5 cm in 
diameter, and developed a large necrotic core. More ne¬ 


crosis was present in the T47D 
tumors than in MCF7 tumors 
having similar age and size. 

Typical image signal-to- 
noise for viable tumor ranged 
from 60:1 on b = 0 s/mm^ im¬ 
ages to 40:1 on b = 754 s/mm^ 
images (where the noise was 
estimated from the standard 
deviation in a ROI in the back¬ 
ground). Areas of necrosis 
could be identified as having 
increased signal intensity on 
the r 2 -weighted images (15, 
16), and reduced signal inten¬ 
sity on the higher diffusion- 
weighted images. MCF7 and 
T47D tumors had similar ap¬ 
pearances for all imaging se¬ 
quences. Representative im¬ 
ages are shown in Fig. 2. 
Figures 2a and 2b show a 
4-week-oid MCF7 tumor (M3), 
where Fig. 2a is a high-resolu¬ 
tion Tg-weighted image, and 
Fig. 2b is a lower resolution 
diffusion-weighted image (b = 
636 s/mm^). Figures 2c and 2d 
show a 4-week-old T47D tu¬ 
mor (T5) for comparison. 

Motion artifacts due to 
breathing were visible in the 
upper abdominal region, but 
very little phase-smearing was 
observed in the region of the 
tumor. Motion of the tumor oc¬ 
curring between scans was 
limited to a few pixels over the 
entire set of b values. This mo¬ 
tion was qualitatively evalu¬ 
ated by outlining the tumor on 
the b = 0 image, and displaying this outline on subse¬ 
quent images. 

Mean ADCs were calculated for viable and necrotic 
regions within the tumors to determine the specificity of 
the ADC for characterizing the tissue. Figure 3a shows 
typical ROIs defined on a high-resolution T 2 -weighted 
image in water, viable tumor, and necrotic tumor for a 
3-week-old T47D tumor (T 2 ). A histological section 
stained with H&E corresponding to this image is shown 
in Fig. 3b. The beginning of necrotic regions near the 
center of the tumor is evident even though this tumor is 
relatively young. A semilog plot of the mean signal in¬ 
tensities within these ROIs is shown in Fig. 3c. 200 X 
magnification views of hematoxylin-eosin-stained histo¬ 
logical sections corresponding to these regions are shown 
in Figs. 3d and 3e. Figure 3d shows a region of viable 
tumor and Fig. 3e corresponds to the necrotic center of 
the tumor. 

ROIs were defined in viable and necrotic tissue on as 
many slices as possible (as selected on the T 2 -weighted 
images) and mean ADCs for these regions are shown in 
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FIG. 3. (a) Ts-welghted image of a 3-week-old T47D tumor (72). Typical ROIs are shown outlining 
regions containing water, viable tumor (white), and necrosis (black), (b) Low-magnification view of 
a H&E-stained histological section corresponding to this image. Areas of necrosis are marked with 
arrows, (c) Mean signal intensities within these ROIs for b = 0 to 754 s/mm^. (d)-(e) 200 x 
magnification views of the H&E-stained sections corresponding to viable tumor and necrosis, 
respectively. 


Table 1. The standard error in the mean for these mea¬ 
surements is reported as the experimental uncertainty. 
Water ADCs were calculated as a reference for each tu¬ 
mor, and reflect the range of temperatures within the 
bore for the experiments. The signal decay from the water 
standards was highly linear. No attempt was made to 
correct measured ADCs in tumors for temperature, since 
the tumor is assumed to be affected little by the ambient 
temperature. The results show a well-defined separation 


between ADCs of necrosis and 
viable tissue for all T47D tu¬ 
mors [P < 0.02 using a Mann- 
Whitney Rank Sum Test). Since 
easily identifiable necrotic re¬ 
gions were present only on two 
MCF7 tumors, no statistical test 
was performed to compare be¬ 
tween necrotic and viable tu¬ 
mor for MCF7. There was no 
statistically significant differ¬ 
ence in calculated ADCs of via¬ 
ble txunor between MCF7 tu¬ 
mors and T47D tumors [P > 
0.05 using a standard t test). 

ADC maps were calculated 
from a straight-line fit of the nat¬ 
ural logarithm of the signal in¬ 
tensity to the diffusion-weight¬ 
ing b factor on a pixel-by-pixel 
basis. These maps show signifi¬ 
cant contrast between areas of 
necrosis and areas of viable 
tumor, and provide a highly 
detailed representation of gross 
tumor morphology. Representa¬ 
tive ADC maps are shown in 
Fig. 4. Figures 4a and 4b show 
one slice firam an 8-week-old 
MCF7 tumor (M4), which has 
begun to develop small regions 
of necrosis near its center. Fig¬ 
ure 4a is an ADC map, while 
Fig. 4b shows a histogram with 
the distribution of ADCs within 
the tumor. The results from an 
older, 11-week-old MCF7 tumor 
(M5) with a large necrotic cen¬ 
ter are shown in Figs. 4cand 4d. 
Two distinct peaks correspond¬ 
ing to the viable and necrotic 
areas are visible in the histo¬ 
gram for this tumor. These ADC 
maps show good correlation 
with gross morphological fea¬ 
tures of the tumors, as observed 
on the histological slices. No 
consistent biexponential signal 
decay could be observed within 
the tumor, even in regions 
showing high vascularity on 
histology where the contribu¬ 
tion to the signal from blood 
flowing in the microcirculation would be expected to be 
strongest. In many pixels, a signal increase was observed at 
low b values, likely indicating that small macroscopic mo¬ 
tions had occurred between scans. 

DISCUSSION 

Tumors from both cell lines grew faster and larger in the 
mammary fat pad than had previously been observed in 
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Table 1 

Measured Apparent Diffusion Coefficients for Regions of Interest in T47D and MCF7 Human 
Breast Tumors 


Apparent Diffusion Coefficients 


Tumor 

uiameier 

(cm) 

Water (10“^ 
mm*/s) 

Necrosis 
(10"^ mm Vs) 

Viable tumor 
(10“^ mm Vs) 

T47D 

T1 (3 weeks)® 

1.0 

2.10 ± 0.07*’ 

2.00 ± 0.07 

0.68 ± 0.03 

T2 (3 weeks) 

1.0 

2.55 ± 0.02 

1.46 ± 0.06 

0.69 ± 0.06 

T3 (4 weeks) 

0.4 

2.52 ± 0.04 

1.57 ± 0.08 

0.49 ± 0.01 

T4 (4 weeks) 

0.5 

2.35 ± 0.02 


0.51 ± 0.07 

T5 (5 weeks) 
Mean value 

1.0 

2.36 ± 0.03 

1.49 ± 0.12 
1.63 ± 0.12 

0.49 ± 0.04 
0.57 ± 0.05 

MCF7 

M1 (4 weeks) 

0.7 

2.34 ± 0.18 


0.47 ± 0.02 

M2 (4 weeks) 

0.8 

2.49 ± 0.07 

1.26 ± 0.08 

0.44 ± 0.05 

M3 (4 weeks) 

1.0 

2.32 ±0.10 


0.48 ± 0.04 

M4 (8 weeks) 

1.5 

2.36 ± 0.08 


0.49 ± 0.02 

M5 (11 weeks) 
Mean value 

2.5 

2.25 ± 0.05 

1.30 ± 0.12 

1.28 ±0.02 

0.43 ± 0.04 
0.46 ± 0.01 


® tumor age refers to the time from day of implantation to the day of measurement. 

The standard deviation of pixel intensities within the RO! is reported as the measured uncertainty For all ROIs 
n > 25 pixels. ’ 
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FIG. 4. (a) ADC map of an 8-week-old MCF7 tumor (M4), which has begun to develop small 
regions of necrosis near its center (areas of high signal intensity), (b) Histogram showing the 
distribution of ADCs within the tumor. (c)-(d) The results from an older, 11-week-old MCF7 tumor 
(M5) with a large necrotic center, (c) ADC map, (d) ADC histogram. 


the flank. MCF7 tumors grew 
faster and to a larger size than 
corresponding T47D tumors in 
this study. We hypothesize 
that this is due to better nutri¬ 
ent and oxygen delivery, as a 
result of increased vascularity 
in MCF7 tumors. More blood 
vessels were apparent on his¬ 
tology in MCF7 tumors than 
T47D. In general, more necro¬ 
sis was present in T47D tu¬ 
mors, and even highly necrotic 
regions in the MCF7 tumors 
contained more cellular de¬ 
bris, and islands of viable tu¬ 
mor cells, indicating more re¬ 
cently formed necrosis. 

Measured ADCs were highly 
specific for necrotic and viable 
tumor in both MCF7 and T47D 
tumors. The higher diffusion 
rates in necrotic areas reflect 
increased volume of extracel¬ 
lular space and increased mo¬ 
bility of water molecules in 
the extracellular compart¬ 
ment. The tortuosity of the ex¬ 
tracellular compartment is re¬ 
duced in necrosis compared 
with viable neoplastic tissue 
due to the decreased density 
of cells. Additionally, the rel¬ 
ative volume fraction of water 
in the extracellular space is 
greatly increased in necrotic 
regions due to cell lysis. In ad¬ 
dition to being specific for 
identification of necrotic re¬ 
gions, it is plausible that the 
measured ADC is an indica¬ 
tion of the stage of the necro¬ 
sis, In initial stages of the for¬ 
mation of necrosis, the cells 
and cytoplasmic compart¬ 
ments swell and the nucleus 
may shrink and become frag¬ 
mented. After disruption of 
the cell membrane, the cell 
contents are released and cel¬ 
lular details are lost. This re¬ 
sults in the appearance of 
sparsely distributed nuclear 
dust, membrane fragments 
and proteins. In addition, in 
most cases we could observe 
islands of still viable cells i 
within the necrotic regions. 
However, the spatial variation 
of these morphologically dis¬ 
tinct features occurs on a scale 
much smaller than the slice 
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thickness of 1 mm. Therefore, the measured ADC for a 
pixel may reflect water mobility at all stages of necrosis, 

L well as in viable tumor. 

Although no statistically significant difference was ob¬ 
served between viable tumor in MCF7 and T47D tumors, 
jhe measured ADCs in necrotic regions in T47D tumors 
were higher than in corresponding necrotic regions in 
MCF7 tumors. This is consistent with what is expected 
from partial volume effects since the necrosis found in 
T47D tumors was, in general, more progressed than that 
found in MCF7 tumors and is therefore expected to be 
ynore edematous. 

ADC maps showed relatively uniform signal intensi¬ 
ties within areas of necrotic and viable tumor, and excel¬ 
lent contrast was observed between necrotic regions 
within the tumors and areas of viable neoplastic tissue. 
The maps compared well to gross tumor morphology as 
evidenced by histological sections corresponding to the 
images. A more detailed comparison to histology could 
not be made because of the difference of a factor of 200 in 
slice thickness between the MR images and histological 
slices. The existence of small macroscopic motions of the 
tumor occuring between scans during the long examina¬ 
tion times precluded an investigation of IVIM effects for 
this data. 

The results suggest that diffusion-weighting imaging, 
with an appropriate choice of weighting parameters 
could provide a valuable noninvasive technique for char¬ 
acterizing tumor response to therapy. Currently, locally 
advanced cancer is managed by using a combined mo¬ 
dality treatment including multiagent chemotherapy, 
surgery, and radiation. Unfortunately, the response of the 
tumors to these therapies is highly variable, and not 
well-predicted, making individual treatment planning a 
very difficult problem. Tumor response to therapy is 
evaluated conventionally by a combination of clinical 
examination and x-ray techniques. However, tumor re¬ 
sponse in the form of cell death and necrosis precedes 
changes that can be detected by using these methods, and 
in many cases, a large degree of necrosis can occur with 
no detectable change in size or radiographic appearance. 
With the development of specialized gradient coil hmd- 
ware to provide adequate diffusion weighting (18), it is 
likely that diffusion imaging by using MR can be used in 
a clinical setting to more efficiently monitor tumor re¬ 
sponse to therapy. 

CONCLUSIONS 

We have applied MR diffusion imaging techniques to hu¬ 
man breast cancer implanted in nude mice. The results of 
this study indicate that the ADC as measured by using a 
PGSE sequence is specific for necrosis and viable tumor. 
Maps of the ADC reflected well the gross morpholo^cal 
features of the tumors as revealed by histological examina¬ 
tion. The results indicate potential for characterizing solid 
tumors by using diffusion imaging, and may have implica¬ 
tions for monitoring tumor response to therapy. 
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Practical Design of a High-Strength Breast Gradient Coil 

Cynthia F. Maier, Hristo N. Nikolov. Kenneth C. Chu, Blaine A. Chronik. Brian K. Rutt 


A high-strength three-axis locai gradient coii set was con¬ 
structed for MRi of the breast. Gradient fieids with good uni¬ 
formity (<10% deviation from the desired gradient) over most 
of the voiume required for breast imaging were generated 
vvith efficiencies of up to 3.3 mT/m/A. The coiis will allow 
diffusion breast imaging in clinically acceptable examination 
times. The electrical design, water cooling system, and fabri¬ 
cation techniques are described. Preliminary tests of the coil 
included images of a grid phantom and diffusion measure¬ 
ments in a short-T 2 agarose gel phantom. 

Key words: breast imaging; diffusion imaging; gradient coiis. 
INTRODUCTION 

In this paper, we present a specialized gradient coil design 
for breast imaging and demonstrate the utility of this coil for 
performing MR diffusion imaging with short-duration, 
high-strength gradient pulses. Sensitization of an MR pulse 
sequence to motion can be achieved using conventionally 
available gradient hardware; however, measurement of wa¬ 
ter diffusion (or slow flow in microvessels) with this hard¬ 
ware requires long echo times to achieve large gradient 
moments and therefore suffers from significant signal loss 
due to T 2 relaxation and system nonidealities. In addition, 
increasing the time between motion-sensitizing gradient 
pulses allows more macroscopic motion to occur in this 
interval, resulting in pixel misregistration, motion-induced 
ghost artifacts, and difiusion/flow measurement artifacts. 
Higher-strength gradients will allow adequate sensitization 
to small motions (i.e., large gradient moments) with very 
short, high-amplitude gradient pulses and minimal time 
delay between these pulses. 

METHODS 

The gradient coil assembly was designed in a cylindrical 
geometry, with each gradient coil built onto a separate 
cylindrical former. The three-axis coil set was composed 
of concentrically nested coils, oriented vertically such 
that their common axis was transverse to the main field 
as shown in Fig. l. Space restrictions within the magnet 
bore require that the coils be as short as possible in the 
vertical direction to allow a patient to lie comfortably 
above the top edge of the cylinders. Unfortunately, de¬ 
creasing the height of the coils (relative to their diame- 
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ters) makes the creation of uniform gradients within the 
coils much more difficult. To satisfy the competing de¬ 
mands of good gradient uniformity and height restric¬ 
tions, we used an aspect ratio of approximately 1:1 for 
our cylindrical coils. The height of the complete coil 
assembly was 16 cm, leaving 34-cm clearance in the 
vertical direction in our General Electric (GE) Signa Ho¬ 
rizon scanner (General Electric Medical Systems, Mil¬ 
waukee, WI), and the largest coil diameter was 21.6 cm. 

To generate the winding densities required for high gra¬ 
dient strengths, relatively small-diameter wire (AWG no. 
20) vvas used to construct the coils. The resistivity of this 
wire is high, requiring an active cooling system to remove 
heat and maintain temperatures inside the coils within a 
range comfortable for the patient. Water-cooling jackets 
with very thin walls were sandwiched between the gradient 
coils, so that resistive heat could be removed efficiently 
while maintaining the compact construction of the coil 
^sembly. A dedicated RF coil and shield were mounted 
mside the innermost gradient coil with an approximately 
1.3-cm space between the RF coil and closest gradient coil. 
The entire assembly is depicted in cross-section in Fig. 2. 

Electrical Design 

The highly efficient quadrupolar cvurent distribution (1—3) 
vvas used to provide gradients in along both the x and z 
directions (G^ and GJ. This current density varies azimuth- 
ally as cos(20) and is entirely directed along the axial di¬ 
rection on the surface of the cylindrical former. Discrete 
wire positions approximating the continuous quadrupolar 
current density were calculated by sampling the current 
density cos(20) at e,- such that cos(2e^+i) “ cos(2e,) was 
constant. A schematic diagram showing this ideal wire 
distribution is given in Fig. 3a. The current pattern in Fig. 
3a is nonphysical in that the current is discontinuous at the 
ends of the coil, and therefore, connecting paths between 
the wires (“return wires”) must be provided. In most imag¬ 
ing applications in which local gradient coils have been 
used previously, the volume to be imaged has been located 
at the center of the coil, and placing return paths near the 
ends of the cylinder had a minimal effect on the desired 
field within the object. For breast imaging applications, the 
required location for the region of good gradient uniformity 
is near the top edge of the coils and a conventional return 
wte arrangement would corrupt the gradient field unifor¬ 
mity within this volume. We chose to place the return wire 
components instead on horizontal plates at the tops and 
bottoms of the coils. Qose to the cylinders on the plates, the 
return paths were extended along the z direction for at least 
5 cm (no is generated by ciurent in the z direction). The 
bottom set of return wires for the coil were strung across 

the bottom of the cylinder in the z direction. The quadru- 
pole current distribution for G, is rotated 45“ with respect 
to G^, and the return wires for G^ could not be strung in the 
z direction; therefore, these wires were located on a hori¬ 
zontal plate. 
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FiG. 1. Schematic showing orientation of local breast gradient 
coll in magnet. 

The design of the Gy gradient makes use of the fact that 
V X = 0 throughout any region in which the current 
density is 0. This means that a transverse gradient coil 
with its axis along z (the conventional orientation) can be 
turned 90® about the x axis so that its axis is along the y 
direction and used to produce a gradient in along the 
y direction. We designed a transverse gradient coil in the 
conventional orientation using the target field method (4, 
5) and then rotated the coil to produce Gy Although a 
wide range of breast sizes exist, we defined an average 
volume of interest (VOI) for breast imaging (by compari¬ 
son with commercially available breast RF coils) as a 
cylindrical volume with length equal to 60% of the coil 
height (11.5 cm), a diameter equal to 90% of the RF coil 
diameter (13.5 cm), and location flush with the top sur¬ 
face of the gradient coil assembly. The target field coil 
was designed to have its region of good gradient unifor¬ 
mity centered on this VOL This design consisted of only 
two sets of closed current loops (Fig. 3b). For the target- 
field coil in its conventional orientation with its axis 
along z, the desired gradient is created almost entirely by 
the current in half of the wire loops only (the “primary” 
component of the loop as indicated in Fig. 3b), and return 
wires are an intrinsic part of the design (labeled “return” 
in Fig, 3b). For this coil, only one set of “return paths” is 
required for current continuity, and in the vertical orien- 



10.82 cm 


FIG. 2. Cross-section of complete gradient coll assembly show¬ 
ing location of RF coll and water-cooling modules. 



a 



b 


Primary 

Wires 


Return 

Wires 


FIG. 3. (a) Quadrupolar current density with 56 windings, (b) Tar- 
get-fleld-design transverse coil. The “primary wire” component is 
almost entirely responsible for generation of the transverse gradi¬ 
ent for the conventional coll orientation with its axis along z. 


tation, these wires are located at the bottom of the cylin¬ 
der. In our coil design, such a large number of loops was 
required to generate the desired gradient efficiency that 
these return wires would have occupied a large fraction 
of the space available on the surface of the cylinder. Due 
to the strict limitations on coil height, we believed that 
this space could be used more efficiently if only the 
primary component of the coil design was present on the 
cylinder itself and the return paths were located on a 
horizontal plate as for the quadrupole coils. A complica¬ 
tion that arises from truncating the wire pattern in the 
vertical dimension in this way is that the reciprocity 
relations from Maxwell’s equations no longer hold for a 
nonphysical cmrent distribution (i.e., V x B 0 for a 
current distribution with nonzero divergence). This 
means that the gradient field calculated for the truncated 
design is not identical in the conventional and transverse 
orientations unless the return wires are included in the 
calculation. We shall see in the Results section that the 
truncated coil in the transverse orientation has both an 
increased gradient efficiency compared to the full design 
for the same coil height and, additionally, the region of 
acceptable gradient uniformity (<20% deviation from 
uniformity) is significantly larger for the truncated coil. 

One undesirable side effect that arises from locating 
the return wires on horizontal surfaces is that the torque 
on the coils from interaction with the main static field is 
greatly increased. This effect is significant. For our 1.5-T 
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magnet, current is only supplied when the coils are po¬ 
sitioned near the central, highly homogeneous region of 
the magnet (typical clinical imaging systems have a 50- 
cm-diameter spherical volume of better than 10 ppm Bq 
uniformity) and the upward (downward) force on the +z 
(-z) half of the horizontal plate for the coil is approx¬ 
imately 15 N/A (3.5 Ib/A) increasing linearly with ap¬ 
plied current. As the gradients are switched, the effect of 
this force is to create an alternating torque of approxi¬ 
mately 3 Nm/A in size. For many pulse sequences, the 
direction of the current is switched rapidly, resulting in 
a torque that alternates direction rapidly and causes sig¬ 
nificant vibrational stress on the coils. For our prelimi¬ 
nary experiments, a temporary coil support was con¬ 
structed from wood- and glass-reinforced phenolic. The 
horizontal plates were rigidly attached to the support, 
and the support structure itself served to mount the coil 
assembly inside the bore of the imager, thereby prevent¬ 
ing any bulk coil motions that could arise from the mag¬ 
netic forces on the wires. The wires themselves were 
potted securely in the milled grooves using epoxy and 
were prevented from vibrating loose from these grooves 
by the next adjacent layer of cooling channels, since all 
layers were tightly sandwiched together. 

Coil Construction 

To produce high gradient efficiencies, it is advantageous 
to use the maximum number of wires possible. This 
number is constrained in an obvious way by the wire 
diameter, minimum spacing between wires, and circum¬ 
ference of the cylindrical former. To obtain the desired 
gradient efficiencies for these coils, it was necessary to 
use relatively small wires (AWG no. 20) with a diameter 
of —1 mm and an interwire spacing of 0.5 mm. The 
(outermost) coil was wound with a single length of wire 
such that 216 axially directed segments generated the 
required current density, and for the (innermost) coil, 
200 axially directed segments were wound. For the mid¬ 
dle coil (Gy), constraining the closest separation of the 
wires to be 0.5 mm allowed 48 loops per side. For ail 
three coils, the wire paths were milled (numerically con¬ 
trolled) into the surface of polycarbonate/acrylic formers. 
Wire was then wound into these grooves and epoxied in 
place. Figure 4a shows the constructed G^ coil, complete 
with horizontal plates after winding with (dielectric) in¬ 
sulated copper wire. The upper plates were curved 
downward at the superior and inferior ends to provide a 
more comfortable structure for patient support. Figure 4b 
shows the middle coil (G^ before assembly of the com¬ 
plete gradient coil unit. 

The resistivity of the small-diameter wire used to con¬ 
struct the coils is relatively high, increasing from 33.3 
mft/m at room temperature to 40.5 mfl/m at 75®C (our 
specified upper limit on temperature at the wire). The 
additional length of wire required to extend the return 
paths away from the main cylinders also contributed to 
high coil total resistances. Since high root-mean-squared 
(RMS) currents are required for the desired imaging ap¬ 
plications, the power dissipated by the coils is significant 
{-^1250 W for one coil @ 25A RMS). In a mock-up of one 
quadrant of a coil, when a 25A RMS current was applied 


with no cooling, the temperature of the quadrant rose 
from 10®C to 100°C in 2.5 minutes. With air cooling only, 
a temperature rise from 10°C to 35°C in 4 minutes was 
observed. A more efficient cooling system was therefore 
required to ensure patient safety. 

Water is one of the most efficient and economical 
cooling fluids available. Its heat capacity (C^ = 4.18 
J/g/®C) is approximately two times higher than other 
commonly available cooling fluids such as hydrocarbon 
oils (e.g., ethylene glycol). Assuming that all heat gener¬ 
ated is transferred to the cooling water and the maximum 
allowable temperature rise of the water is 5®C, a water 
flow of 3.6 liter/min is required to remove 1250 W of 
heat. Conventionally, water-cooling channels for gradi¬ 
ent coils are constructed from copper tubing. For a com¬ 
pact gradient coil set, the competing demands of space 
limitations and the requirement of a large cooling chan¬ 
nel cross-section to ensure adequate volume flow rates 
necessitated the design of a novel cooling system. We 
designed this cooling system in separate modules, each 
module corresponding to a quadrant of a single gradient 
coil. Each module consisted of a skeleton of cooling 
channel spacers that were milled from a 3.0-mm-thick 
polycarbonate sheet. This was covered on both sides 
with a very thin (0.25-mm-thick) polycarbonate film (GE 
Plastics, General Electric Company, Pittsfield, MA) to 
allow close thermal contact of the wires with the cooling 
water while ensuring electrical insulation. The film was 
bonded to the skeleton using ultraviolet-curing polyvinyl 
chloride-based bonding adhesive (Loctite 3311, Loctite 
Corporation, Mississauga, Canada). The tensile strength 
(to break) of the polycarbonate film was rated at 9000 psi, 
according to American Standards of Testing Materials 
(ASTM) D882, and the burst strength was 200 psi (ASTM 
D774). The modules were tested for flow uniformity and 
resistance to leaks by flowing water at 13 liter/min for >5 
hours before permanent incorporation into the gradient 
coil structure. After incorporation into the gradient coil 
set, the modules were connected in parallel to lower the 
fluid resistance of the network of cooling channels and 
thereby allow increased water flow rates. 

This design allowed very close thermal contact be¬ 
tween the wires and water in the cooling channels (while 
maintaining electrical isolation) and also provided a 
large surface area over which heat could be transferred to 
the water. The ability to test the cooling channel integrity 
separate from the gradient coil structure was vital in 
ensuring no electrical shorting due to water leakage. An 
additional benefit of constructing the cooling system 
from polycarbonate is that no electrical currents can be 
created in this material, and it is therefore not susceptible 
to eddy current generation like some copper systems. 
Since the thickness' of the cooling channel module was 
only 3.5 mm, the maximum total wall thickness of the 
gradient coil structure was less than 3 cm. The assembled 
gradient coil structure, including RF coil and water con¬ 
nections, is shown in Fig. 5. 

RF Coil 

Several cylindrical RF coil designs exist that can be used 
in a transverse geometry for imaging the breast (i.e., with 
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FIG. 4. (a) Gjf quadrupole coil with 216 windings, (b) coil showing placement of the “return wires” on a horizontal plate, 

the coil axis perpendicular to the main magnetic field proximity of the RF coil to the gradient coils results m 

Bq). These coils, when oriented vertically, produce very strong inductive coupling and detoing of the RF 

fields directed along either the horizontal x direction coil (given that the gradient coil radius is constrained to 

(e.g., saddle coils (6)1 or along the vertical direction (e.g., be as small as possible for maximum efficiency). To 
solenoidal coils (7, 8)). determine the best RF coil design for our gradient coil 

The construction of an RF coil for operation inside our set, two different RF coils were constructed and testecL 

breast gradient coils is a difficult problem, in that the The first coil constructed was a quadrature RF coil with 
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FIG. 5. Assembled gradient coil module showing RF col! and 
cooling channels. 


two separate coils providing the perpendicular com¬ 
ponents. In this arrangement, a simple three-loop sole- 
noidal coil provided the vertical mode, and a saddle coil 
built onto the same former provided the horizontal mode. 
The second coil constructed was a 16-rung high-pass 
birdcage RF coil designed for operation in a linear mode 
(recall that a birdcage coil cannot be used in quadrature 
mode in the vertical orientation because one of the 
modes would be parallel to Bq), We chose this design 
over a saddle-coil arrangement because of the superior 
uniformity of produced by the birdcage coil over the 
volume enclosed by the coil. A split cylindrical RF shield 
was constructed from fine copper screen and located 
immediately inside the smallest diameter gradient coil. 
The coils were tested for coupling to the gradient coils. 
The quality factors, Q, were measured, and the coil per¬ 
formances were compared. 

Diffusion Imaging 

Preliminary diffusion measurements were made using a 
short-T 2 phantom consisting of NiS 04 -doped agarose gel 
(72 = 60 ms at 1.89 T) with a diffusion-sensitized spin- 
echo sequence (pulsed gradient spin-echo sequence (9)). 
Apparent diffusion coefficient (ADC) maps were calcu¬ 
lated from a series of 10 diffusion-weighted images with 
diffusion b factors between 0 and 1000 s/mm^ Each of 
these images was the geometrical average of two images 
with equal diffusion-weighting but with opposite diffu¬ 
sion gradient polarity to remove the effects of diffusion¬ 
weighting cross-terms from interaction between the dif¬ 
fusion and imaging gradients (10). The ADC maps were 
obtained under three conditions: (1) using the standard 
whole-body manufacturer-supplied gradient coils for all 
imaging and diffusion-weighting gradients; (2) using our 
custom coils for all imaging and diffusion-weighting gra¬ 
dients, with the gradient strengths matched to those of 
the whole-body coils by applying an appropriately scaled 
current to the coils; and (3) using our custom coils for all 
imaging and diffusion-weighting gradients, with the im¬ 
aging gradient strengths matched to those of the whole- 
body coils, but taking advantage of the increased gradient 
efficiency to apply higher strength, shorter diffusion¬ 
weighting gradients, thereby shortening the scan TE. In 
the first two experiments, the minimum TE required for 


diffusion weighting of b = 995 s/mm^ was TE = 80 ms 
(with diffusion gradient strength of 22 mT/m). TR was set 
to 600 ms, the matrix size for the scan was 256 X 128, and 
the image size was 8X8 cm, for a total scan time of 2.5 
minutes for each image (5 minutes for each b factor). 
Diffusion gradients were applied along the readout axis, 
X, for all experiments. For the third experiment, the TE 
was decreased to TE = 46 ms by increasing the strength 
of the diffusion-weighting gradient pulses by a factor of 
four (88 mT/m) to obtain much larger b factors with 
shorter diffusion pulses. All other parameters (including 
slice thickness) were matched to the first two experi¬ 
ments. 

RESULTS 

Coil Performance 

The field produced by each of the three coils was calcu¬ 
lated by direct numerical solution of the Biot-Savart Law 
for the wire paths, including return wire contributions. 
Calculated gradient efficiency for all three coils, together 
with the measured resistance and inductance, are sum¬ 
marized in Table 1. Contour plots of calculated gradient 
uniformity for the quadrupole coil are shown in Fig. 6. 
Figure 6a shows the percentage deviation from gradient 
uniformity for a y-x plane through the center of the coils. 
A similar plot is shovm in Fig. 6b for a x-z plane, trans¬ 
verse to the coil axis, at a height equal to the center of the 
VOI (i.e., located 5.9 cm below the top surface of the 
coils). These plots show the excellent gradient unifor¬ 
mity within the VOI that is predicted for the coil. The 
gradient strength decreases near the ends of the coil, 
because the finite length of the coil causes a decrease in 
field strength. Similar results were obtained for the G^ 
gradient coil, with appropriate scaling for the number of 
wires and the coil radius. The performance of the modi¬ 
fied target field Gy coil is shovm in Figs. 6c and 6d. A 
larger deviation from gradient uniformity is present for 
this coil, compared to the and G^ coils. This problem 
is inherent to any design for a Gy coil in the transverse 
orientation, since there is an additional gradient in the y 
direction due to the natural field fall-off that occurs as 
the ends of the coil are approached. To create a Gy coil 
with optimum gradient efficiency, we would like to max¬ 
imize jB^l at both ends of the coil, with B^ a maximum in 
the +z direction at one end of the coil, and a maximum 
in the —z direction at the opposite coil end. The field 
fall-off inherent in this coil orientation works against the 
gradient we desire to create, since it necessarily results in 
a decrease in the magnitude of B^ as the coil ends are 
approached. Instead of working against the natural field 
fall-off toward the ends of the coil, we have designed a 


Table 1 

Coll Performance Characteristics 



Mean 

diameter 

(cm) 

Efficiency 

(mT/m/A) 

Inductance 

(mH) 

Resistance 

(n) 

Gx (quadrupole) 

21.1 

2.9 

754 

2.3 

Gy (fingerprint) 

19.5 

1.7 

829 

2.6 

(quadrupole) 

18.1 

3.3 

1024 

1.9 
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FIG 6 Contour plots showing percentage deviations from gradient uniformity for the coil over (a) an axial plane (x->0 showing the 
entire VOI with the top edge of the plot corresponding to the top surface of the coils. The outer circumference of the VOl Is at 6,75 cm 
as indicated by dashed lines, (b) A coronal plane (x-z) located at the center of the defined VOI, i.e., at ~3 cm above the geometric ^il 
center. In Fig. 6a, the contours are symmetric around y = 0 cm with some small deviations due to return-wire contributions, (c, d) The 
corresponding plots for the constructed Gy coil shown in Fig. 4b. (e, f) The corresponding plots for a conventional transverse coil with 
return wires on the cylinder, as in Fig. 3b. 


coil to produce = 0 at the upper end of the coil, and a 
maximum in toward the bottom end, similar to the 
as 3 mimetric head coil design by Roemer (11)- We used 
the target field approach to optimize the rate at which the 
field falls off toward the upper end of the coil to produce 
a linear variation in Although the region of acceptable 
gradient uniformity for our truncated target field coil is 
not as large as for the quadrupole coils, it compares 
favorably to the performEmce of the full target field de¬ 
sign. We show in Figs. 6e and 6f the corresponding plots 
from a complete 40 loop per side target field coil with the 
same dimensions as our modified coil (i.e., with return 
wire paths at the bottom of the cylinder as in Fig. 3b). The 
region of acceptable gradient uniformity for our trun¬ 
cated target field coil is much larger than for the com¬ 
plete target field coil. In addition, the gradient efficiency 
is better for our truncated coil by a factor of 1.7, due to 


he larger number of windings that could be placed on 
he cylinder for our coil. 

To verify the field simulations, gradient uniformity 
Aras tested by acquiring images of grid phantoms. Devi- 
itions fi'om gradient uniformity are manifest in the im- 
iges as geometric distortions. The phantoms consisted of 
icrylic cylinders filled with square plastic gridding and 
CuS 04 -doped 0.9% saline solution. The grid spacing was 
1.7 X 1.7 cm, whereas the phantoms had an outer diam- 
3 ter of 12.5 cm and a height of 15 cm (extending 3 cm 
higher than the top surface when placed inside the coils). 
The coil performance was tested on a clinical GE Signa 
imager, using our coils for all imaging gradients (with 
gradient strengths matched to the whole-body coils) in a 
standard spin-echo pulse sequence. In general, local gra¬ 
dient coils often require modification of the compensa¬ 
tion network for the gradient coil amplifiers since the 
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coil inductances are quite different from the whole-body 
coils. Because the inductances of our breast coils were 
very similar to the whole-body coils, no modification of 
the network was necessary. Images of a grid phantom are 
shown in Fig. 7 for planes corresponding approximately 
to the theoretical field plots shown in Figs. 6a and 6b and 
are evidence of the good gradient uniformity achieved 
over most of the VOL These images were acquired using 
a spin-echo sequence with TE = 20 ms, TR = 500 ms, 
matrix size = 256 X 256, and 2 NEX. The image size is 
16 X 16 cm, and the VOI vertical extent is marked with 
arrows. 

The first RF coil we constructed was a quadrature coil 
consisting of a solenoidal coil to provide the vertical 
mode and a saddle coil to provide the horizontal coil. 
Unfortunately, the vertical mode was inductively cou¬ 
pled to both the y-gradient coil and to the RF shield and 
could not be tuned to the proton resonant fi-equency at 
1.5 T. We therefore chose to construct a coil that would 
operate in linear mode, with only a horizontal com¬ 
ponent present. We constructed a birdcage coil that cre¬ 
ates a horizontal mode field with a high degree of 
homogeneity over a large volume. The unloaded Q of this 
coil (when located inside the gradient coil structure) was 
80, whereas the Q for loading with a 0.9% saline phan¬ 
tom was 29. 

Although the gradient strengths required for conven¬ 
tional imaging (10 mT/m) can be generated at less than 6 
A using these coils, higher currents are required for more 
demanding applications such as flow imaging and diffu¬ 
sion imaging. To test the effectiveness of the water cool¬ 
ing system, a DC current was applied to the coil with the 
highest total resistance, Gy, and the temperature rise was 
measured for different water flow rates by monitoring a 
thermocouple located on the surface of the wires. The 
temperature was monitored until equilibrium was 
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achieved. Measured temperature rises are shown in Fig. 8 
for DC currents between 0 and 22 A and a cooling water 
flow rate of 8 liter/min. Similar results were obtained for 
other flow rates, at which the equilibrium temperature 
increased with decreasing flow rates but did not change 
significantly from the results shown in Fig. 8 for higher 
flow rates. This behavior has also been demonstrated in 
purely cylindrical gradient coils and was modeled by 
Chu and Rutt (12) using an analysis based on fundamen¬ 
tal heat transfer theory. 

Diffusion Imaging 

ADC maps for a T 2 = 60 ms NiS 04 -doped agar phantom 
are shown in Fig. 9; Fig. 9a shows the results obtained 
using the GE gradient coils, and Fig. 9b shows the results 
obtained using our gradient coils with diffusion gradient 
strengths four times higher than that of the GE coils. The 
mean ADC for all three experiments was the same within 
experimental error: (1.8 ± 0,3) X 10"^ mm^/s for the GE 
coils (±30% peak-to-peak), (1.8 ± 0.3) x 10“^ mm^/s for 
our gradient coils with diffusion gradient strengths 
matched to the GE coils, (1.9 ± 0.2) X 10“*^ mm^/s for our 
gradient coils with diffusion gradient strengths four 
times higher than that of the GE coils and TE = 46 ms. As 
expected, the observed signal-to-noise ratio was in¬ 
creased by nearly a factor of 2 for the TE = 46 ms images 
compared to the corresponding h-factor images at TE - 
80 ms. 


DISCUSSION 

The performance and homogeneity of the gradient coils 
were in excellent agreement with theoretical predictions. 
The volume of acceptable gradient uniformity (<20% 
deviation) extended over most of our defined VOI for all 



FIG. 7. Spin-echo images of a grid phantom acquired (a) for an axial plane (x-y) and (b) for a coronal plane (x-z) with the specialized breast 
quadrupole gradient coils. The input current strength required for the least efficient (Gj,) coll was approximately 14 times less than that 
required for the corresponding GE coil. 
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fig. 8. Temperature response of assembled gradient coils to a 
DC current applied to the coil with the highest resistance (G^) 
monitored at the surface of the G^ coll where the winding density 
was greatest. 


three gradient coils. For this region, the theoretical plots 
of can be used to generate a gradient-warping algo¬ 
rithm for the coils (13) to correct spatial distortions in the 
images. Additionally, the known deviations from gradi¬ 
ent uniformity can be used to correct calculations of 
diffusion-weighting b factors (which will vary spatially 
across the VOI for our coils). Unfortunately, greater de¬ 
viations from uniformity could not be avoided for the top 
1 cm of the VOI near the upper edges of the coils for our 
coil designs. Deviations from uniformity of more than 
20% lead to large image distortions and likely cannot be 
satisfactorily corrected using conventional schemes. 
Many breast lesions occur close to the chest wall; to be 
clinically useful, it is important to image accurately as 
close to the chest wall as possible. Several new gradient 
coils with improved gradient uniformity at the top edge 
of the coil have been reported recently (14, 15). Future 
improvements on our coil design will focus on extending 
the region of usable gradient beyond the top edge of the 
coil (16) as far as possible. For our preliminary clinical 
work with these coils, we will limit our study to patients 
with known lesions occurring in a region of the breast 
that falls within the 20% uniformity region of our gradi¬ 
ent coils. This approach will allow us to implement 
pulse sequences under ideal conditions, (a) to determine 
the best parameters for visualizing breast lesions, (b) to 
determine the minimum required gradient strengths and 
uniformities, and (c) to make an initial assessment of the 
clinical utility of this approach. 

In future coil designs, it may be possible to extend the 
region of usable gradient farther into the cylinder in the 
-y direction by recessing the coils in the patient table or 
by constructing a dedicated breast imaging table, thereby 
allowing higher aspect ratio coils. In this design, the coil 
assembly would be constructed to allow positioning on 


either side of the table, depending on which breast was 
being imaged. A fiberglass or Kevlar sheath covering the 
coil assembly and running the length of the patient table 
would serve the dual functions of isolating the patient 
from the electrical and cooling circuits and providing an 
additional form of mechanical fixation to secure the coils 
on the table. Because the coils are not torque-balanced, it 
will be necessary to anchor the structure firmly, both 
from the point of view of vibration-induced imaging ar¬ 
tifacts and of patient safety. 

The high gradient efficiencies of these coils will allow 
much higher spatial-resolution imaging than can cur¬ 
rently be achieved in a clinical setting, and additionally, 
will provide very large amplitude motion-sensitizing gra¬ 
dients, thus facilitating diffusion imaging (or flow imag¬ 
ing) in small blood vessels with short echo times. T 2 has 
been measured for human breast tissue, with reported 
values ranging between 60 and 150 ms (17,18). Recently, 
we reported T 2 approximately 70 ms for human breast 
tumor (MCF7 human breast carcinoma) growing in a 
mouse mammary fat pad (19). More fibrous tumors can 
have significantly shorter T 2 S. For this reason, initial 
signal loss at b = 0 s/mm^ due to T 2 weighting of the 
image can severely limit the range of b factors for which 
images can be obtained, which has the result of limiting 
the precision of the ADC measurement. Our approach is 
to shorten TE while maintaining diffusion weighting by 
increasing diffusion gradient amplitudes. This gives a 
substantial increase in initial signal intensity for b = 0 
s/mm^, allowing high-quality images to be obtained for 
much higher b factors. Even in cases in which signal loss 
from T 2 weighting is not so severe as to limit ADC pre¬ 
cision, the ability to obtain high diffusion weighting with 
short TE reduces dramatically the T 2 contrast in an im¬ 
age. This will allow acquisition of a diffusion-weighted 
image in which the contrast is determined mainly by 
tissue ADC rather than by a mix of ADC and T 2 . 

A secondary important benefit that can be gained using 
high-efficiency gradient coils is that the ramp time re¬ 
quired to attain imaging gradient strengths (10-25 
mT/m) is decreased from that required by whole-body 
coils by a factor of the efficiency (since the current slew 
rate is similar to whole-body coils, and the required peak 
current is decreased). This will enable echo-planar imag¬ 
ing (and other fast imaging techniques) using conven¬ 
tional gradient amplifiers. Echo-planar imaging tech¬ 
niques have been proven to be much less susceptible to 
motion artifacts than conventional imaging sequences 
(20) and have recently been applied to diffusion imaging 
in vivo (21) with good results. 

For rapid imaging pulse sequences, and for high-am¬ 
plitude bipolar pulses, very rapid changes in gradient 
polarity are required. Although the magnetic field ampli¬ 
tudes drop off quickly in the region above the coils, 
magnetic fields are also generated by the return wires on 
the horizontal plate, and excitation thresholds for cardiac 
stimulation must be considered. As an upper limit for 
one coil, |J?| was calculated for the coil at various 
distances above the return wires at the location of the 
highest wire density. The maximum dl/dt from our gra¬ 
dient coil amplifiers is 200 A in 184 /ms. This produces 
dB/dt =156 T/s at 1 cm above the coil return wire plate 
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FIG. 9. ADC maps of a fg = 60 ms agar phantom obtained (a) using the whole-body coils to provide all diffusion and Imaging gradients, 
TE = 80 ms, and (b) using the breast gradient coils to provide all gradients, where the diffusion gradient strengths were four times the 
whole-body gradient strengths and TE = 46 ms. In all cases, the ADC maps were constructed from 10 Images with b = 0 to 6 = 1000 
s/mm^. 


and dB/dt = 50 T/s at 10 cm, compared to a recom¬ 
mended upper limit of 400 T/s for cardiac stimulation 
(22). The magnitude of these fields is relatively small at 
regions of concern with respect to magnetic stimulation 
such as the heart; however, peripheral nerve stimulation 
has been reported at dB/dt as low as 40 T/s in some 
patients (22). Theoretical field calculations will allow 
prediction of |B| and dB/dt for all pulse sequences before 
implementation, and pulse sequences will be tailored to 
adhere to maximum dB/dt and peak |B| guidelines. Ad¬ 
ditionally, if peripheral nerve stimulation is reported for 
individual patients, dB/dt will be reduced below the 
stimulation threshold by increasing gradient switching 
times. 


SUMMARY 

We have constructed a three-axis gradient coil set for 
imaging human breast tumors. These coils create gradi¬ 
ent efficiencies of up to 3,3 mT/m/A with relatively low 
inductance (^1000 pH), and good gradient uniformity 
over most of the volume required for breast imaging. 
Preliminary diffusion results were obtained in a water 
phantom for b factors as high as ~1000 s/mm^ for TE = 
46 ms. A compact, efficient water-cooling system was 
presented, which maintained coil temperatmes at or be¬ 
low room temperature, even for RMS currents as high as 
20 A. This gradient/cooling module will allow, for the 
first time, clinical breast imaging using diffusion and 
microvascular flow-weighted protocols with high spatial 
resolution. 
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